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Colloids immersed in a critical or near-critical binary liquid mixture and close to a chemically patterned substrate are 
subject to normal and lateral critical Casimir forces of dominating strength. For a single colloid we calculate these 
attractive or repulsive forces and the corresponding critical Casimir potentials within mean-field theory. Within this 
approach we also discuss the quality of the Derjaguin approximation and apply it to Monte Carlo simulation data 
available for the system under study. We find that the range of validity of the Derjaguin approximation is rather large 
and that it fails only for surface structures which are very small compared to the geometric mean of the size of the 
colloid and its distance from the substrate. For certain chemical structures of the substrate the critical Casimir force 
acting on the colloid can change sign as a function of the distance between the particle and the substrate; this provides a 
mechanism for stable levitation at a certain distance which can be strongly tuned by temperature, i.e., with a sensitivity 
of more than 200nm/K. 



o 



I 

C 
O 
o 



> 

(N 

oo 



in 

o 
o 



X 



PACS numbers: 05.70. Jk, 82.70.Dd, 68.35.Rh 



I. INTRODUCTION 

Since the discovery of the Casimir effect in quantum 
electrodynamics^ 2 - it is well-known that the inherent fluctua- 
tions of a medium lead to an effective force acting on its con- 
fining boundaries. In soft matter physics, the analogue of the 
vacuum fluctuations in quantum electrodynamics are the ther- 
mal fluctuations of the order parameter of a fluid. These oc- 
cur on the length scale of the bulk correlation length £, which 
is generically of molecular size. However, upon approach- 
ing a critical point at the temperature T = T c , the correlation 
length 4 increases with an algebraic singularity and attains 
macroscopic values. The confinement of these long-ranged 
fluctuations results in the so-called critical Casimir force act- 
ing on a length scale set by <i; 3 . Since the correlation length di- 
verges as 4 (T — >• T c ) °c \T — T c \~ v , where V is a standard bulk 
critical exponent, the range of the critical Casimir force (and 
therefore its strength at a certain distance) can be controlled 
and tuned by minute temperature changes (see, e.g., Refs. 4 
and 5). The characteristic energy scale of the critical Casimir 
effect is given by ksT c , which allows for a direct measurement 
of the critical Casimir forces, in particular if the critical point 
is located at ambient thermodynamic conditions 6 - 7 . 

The attractive or repulsive character of the critical Casimir 
force can be controlled by suitable treatments of the confining 
surfaces. Generically, the surfaces which confine a binary liq- 
uid mixture preferentially adsorb one of its two components 
(or the gas or liquid phase in the case of a one-component 
fluid). This can be described by effective, symmetry breaking 
surface fields, which lead to a preference for either positive 
[(+)] or negative [(— )] values of the scalar order parameter 
<p, corresponding to the difference between the local concen- 
trations of the two species (or the deviation of the density of 
the one-component fluid from its critical value). The criti- 
cal Casimir force strongly depends on the effective boundary 
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conditions (BC) at the walls (see, e.g., Refs. 181-11 51 and refer- 
ences therein). It is attractive for equal symmetry breaking 
(±,±) BC and repulsive for opposing (±,=f) BC. Inter alia, 
this latter feature qualifies critical Casimir forces to be a tool 
to overcome the problem of "stiction" which occurs in micro- 
and nano-mechanical devices. (The quantum electrodynamic 
Casimir force is typically attractive and thus responsible for 
stiction; turning it to be repulsive requires a careful choice of 
the fluid and of the bulk materials of the confinement 16 .) The 
theoretical description of the critical Casimir forces is partic- 
ularly challenging due to the non-Gaussian character of the 
order parameter fluctuations, which contrasts with the intrin- 
sically Gaussian nature of the low energy fluctuations of the 
electromagnetic field; in addition, the critical Casimir effect is 
also particularly rich as it allows, inter alia, symmetry break- 
ing boundary conditions, which do not occur for electromag- 
netic fields. 

The critical Casimir effect exhibits universality, i.e., the 
critical Casimir force expressed in terms of suitable scaling 
variables depends only on the universality class of the bulk 
critical point and on the type of boundary conditions, whereas 
it is independent of the microscopic structure and of the ma- 
terial properties of the specific fluid medium involved. In our 
present theoretical analysis we focus on the Ising universality 
class which encompasses the experimentally relevant classical 
binary liquid mixtures and simple fluids. 

The existence of the critical Casimir effect has been exper- 
imentally confirmed and its strength has been first measured 
indirectly for wetting flimsier— . The first direct measurement 
of this effect has been performed at the sub-micrometer scale 
for a spherical colloid immersed in a (near) critical binary 
liquid mixture close to a laterally homogeneous and planar 
substrate 6 ' 7 . The corresponding Monte Carlo simulation data 
for the film geometry are in very good quantitative agreement 
with all available experimental data 6,7 ' 21 " — . Theoretical stud- 
ies of the critical Casimir effect acting on colloidal particles 
involve sphericall y 25 " 29 or ellipsoidally 30 shaped colloids ad- 
jacent to homogeneous substrates. 

Besides their wide use as model systems in soft matter 
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FIG. 1 . Sketch of a spherical colloid immersed in a near-critical bi- 
nary liquid mixture (not shown) and close to a (patterned) planar sub- 
strate. The sphere with (b) boundary condition (BC) and radius R is 
located at a surface-to-surface distance D from the substrate and its 
center has a lateral coordinate x = X with the substrate pattern being 
translationally invariant in all other directions. The following four 
different types of substrate surfaces are considered: homogeneous 
substrate rSec. HIII . a chemical step [s; Sec. lIVI . a single chemical 
lane [£; Sec.|V), and a periodically patterned substrate [p; Sec. lVll . 
(Note that for a four-dimensional system, which we also consider, 
this is a three-dimensional cut of the system, which is invariant along 
the fourth direction; the sphere thus corresponds to a hypercylinder 
in four dimensions.) For later reference, the box on the left side 
summarizes the definitions of the various scaling variables which the 
scaling functions of the critical Casimir force depend on for the listed 
geometrical configurations. On the right, (a), (a>), (ai), (tfi), and 
(02) indicate the boundary conditions corresponding to the various 
chemical patterns. 



physics, colloids have applications at the micro- and nanome- 
ter scale. In this context, they are widely used in micro- and 
nano-mechanical devices. Therefore, one may utilize the crit- 
ical Casimir forces acting on colloids because their strength 
and their direction can be tuned in a controlled way. Suitably 
designed chemically or geometrically structured substrates 
generate lateral critical Casimir forces acting on colloidal 
particles 3 -^ - —. Current techniques allow one to endow solid 
surfaces with precise structures on the nano- and micrometer- 
scale. Hence, the critical Casimir effect can be used to create 
laterally confining potentials for a single colloid, which can 
be tuned by temperature 32 . 

Recently, the critical Casimir potential of a colloid close to 
a substrate with a pattern of parallel chemical stripes with lat- 
erally alternating adsorption preference has been measured 32 . 
In our corresponding theoretical study 31 , we have calculated 
the normal and lateral critical Casimir forces acting on a col- 
loid close to such a patterned substrate as well as the corre- 
sponding potentials. We have used our theoretical predictions 
for the universal scaling functions of the critical Casimir po- 
tential in order to interpret the available experimental data in 
Ref. [3J- It has turned out that an agreement between theory 
and experiment can be achieved only if one takes into account 
the geometrical details of the chemical substrate pattern. This 



demonstrates that the critical Casimir effect is very sensitive 
to the details of the imprinted structures and that it can resolve 
them. 

Here we generalize our previous analysis 3 -! to various sub- 
strate patterns. In particular we study the critical Casimir ef- 
fect for a f/zree-dimensional sphere close to a homogeneous 
substrate [Sec. HITll . a chemical step [Sec lIVL a single chemi- 
cal lane [Sec.[V], and periodic patterns of chemical stripes of 
alternating adsorption preference [Sec. IVIII [see Fig. [Q. For 
completeness, we also consider a cylinder which is aligned 
with the chemical pattern [Sec. lVHl . We provide quantitative 
predictions for the scaling functions of the critical Casimir 
forces, pursuing a two-pronged approach: (i) We calculate 
the force using the full three-dimensional numerical analysis 
of the appropriate mean-field theory (MFT). (ii) We use the 
so-called Derjaguin approximation (DA) based on the scal- 
ing functions for the critical Casimir force in the film geom- 
etry either obtained analytically within MFT— or obtained 
from Monte Carlo simulations 2 ^ 3 -, which allows us to pre- 
dict the critical Casimir force in the physically relevant three- 
dimensional case. Inter alia, we determine the range of valid- 
ity of the DA within MFT, which provides guidance concern- 
ing its applicability in three spatial dimensions d = 3. This is 
an important information because presently available Monte 
Carlo simulations are far from being able to capture complex 
geometries 2 ^ 3 .. 

Currently, the possibility of realizing stable levitation of 
particles by means of the electrodynamic Casimir forces has 
been the subject of intense theoretical investigation 36-4 !. Our 
results presented in Sees. [VI] and |VII| show that for suitable 
choices of the geometry of the chemical pattern of the sub- 
strate, the critical Casimir forces can be used to levitate a col- 
loid above the substrate at a height which can be tuned by 
temperature. This levitation is stable against perturbations be- 
cause it corresponds to a minimum of the potential of the crit- 
ical Casimir force acting on the colloid. 

In Sec. [II] we briefly introduce the necessary terminology 
related to finite-size scaling and we discuss briefly the corre- 
sponding MFT. Section|III]is devoted to the well-studied case 
of a colloid close to a homogeneous substrate. (In d = 4, as ap- 
propriate for MFT, the three-dimensional colloid is extended 
to the fourth dimension as a hypercylinder, for which we also 
present the results of our analysis.) As mentioned above, the 
various patterns and setups are considered in Sees. IIVT4VHI 
We conclude and summarize our findings in Sec. IVIIII Cer- 
tain important technical details concerning the calculation of 
the Derjaguin approximation are presented in the Appendices 

EHd] 



II. THEORETICAL BACKGROUND 

A. Finite-size scaling 

According to the theory of finite-size scaling, the normal 
and lateral critical Casimir forces and the corresponding po- 
tentials can be described by universal scaling functions, which 
are independent of the molecular details of the system but 
depend only on the gross features of the system, i.e., on the 
bulk universality class (see, e.g., Refs.[8|and[9|and references 
therein) of the associated critical point. Here, we focus on the 



Ising universality class (which is characterized by a scalar or- 
der parameter 0) in spatial dimensions d — 3 and d — 4. In ad- 
dition, the critical Casimir force depends on the type of effec- 
tive boundary conditions at the walls, which we denote by (a) 
and (b), and by the geometry of the confining surfaces 4 - 2 - - ^ 4 -. 
Note that (a) and (b) can represent the various symmetry pre- 
serving fixed-point BC (the so-called ordinary, special, peri- 
odic, or antiperiodic boundary conditions 8 - 9 ) in addition to the 
symmetry breaking cases (±) we are mainly interested in, and 
which describe the adsorption of fluids at the confining walls. 
Inspired by the experiments described in Ref. 32 we con- 
sider binary liquid mixtures with their consolute critical point 
approached by varying the temperature T towards T c at fixed 
pressure and critical composition. We first study the film ge- 
ometry in which the fluid undergoing the continuous phase 
transition is confined between two parallel, infinitely extended 
walls at distance L. According to renormalization group the- 
ory the normal critical Casimir force f( a M P er un it area which 
is acting on the walls scales as 10 



f(a,b)(L,T)=k B Tjjk {aJ , ) (sign(t)L/^ 
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where (a,b) denotes the pair of boundary conditions (a) and 
(b) characterizing the two walls. The scaling function kr a b) 
depends only on a single scaling variable given by the sign 
of the reduced temperature distance t from the critical point 
(± for t ^ 0) and the film thickness L in units of the bulk 
correlation length E,±(t -> ± ) = t,^\t |~ v , where v ~ 0.63 in 
d = 3 and V = 1 /2 in d = 4 45 . (Clearly, one has f (ah) (L,T) = 
f(b,a)(L,T).) Positive values of t, t > 0, correspond to the 
disordered (homogeneous) phase of the fluid, whereas neg- 
ative values of t, t < 0, correspond to the ordered (inhomo- 
geneous) phase, where phase separation occurs. Typically, 
the homogeneous phase is found at high temperatures, and 
one has t = (T — T c )/T c , However, many experimentally rel- 
evant binary liquid mixtures exhibit a lower critical point, 
for which the homogeneous phase corresponds to the low- 
temperature phase and one has t = — (T — T C )/T$L. The 
two non-universal amplitudes £f of the correlation length 
are of molecular size and characterized by the universal ratio 
^/^-^l.9md^3^^md^+/^- = V2md^4 41 ^ ± is 
determined by the exponential spatial decay of the two-point 
correlation function of the order parameter (j) in the bulk. 

At the critical point T = T c , the correlation length diverges, 
£± — s- oo, and the scaling function of the critical Casimir force 
acting on the two planar walls attains a universal constant 
value referred to as the critical Casimir amplitude 8 - 9 : 
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Away from criticality, the critical Casimir force decays ex- 
ponentially as a function of L/£±. For the specific case of 
symmetry breaking BC a,b e {+, — } and for t > one ex- 
pects for L/S.+ >la pure exponential decay of ft ±+ \ (see, 
e.g., Refs. 0, [H and SI and footnote 3 in Ref. 31, i.e., a 
decay without an algebraic prefactor to the exponential and 
without a numerical prefactor to L/|+ in the argument of the 
exponential) corresponding to 
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where A± are universal constants 7 . Note that, in the absence 
of symmetry-breaking fields inside the film, the scaling func- 
tions for (+, +) BC are the same as for (— ,— ) BC. 



B. Mean-field theory 

The standard Landau-Ginzburg-Wilson fixed-point effec- 
tive Hamiltonian describing critical phenomena of the Ising 
universality class is given by 4 ^ 4 - 3 - 
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where (j) (r) is the order parameter describing the fluid, which 
completely fills the volume V in ^-dimensional space. The 
first term in the integral in Eq. dU penalizes local fluctuations 
of the order parameter. The parameter T in Eq. (0]i is propor- 
tional to t , and the coupling constant u is positive and provides 
stability of the Hamiltonian for t < 0. The mean-field order 
parameter profile m := u l / 2 (<p) minimizes the Hamiltonian, 
i.e., 8J$?[<p]/8<p\ ( , =u -i/2 m = 0. In the bulk the mean-field or- 
der parameter is spatially constant and attains the values (<j>) = 
±a\t\P for t < and (0) = for t > 0, where, besides |q" , a is 
the only additional independent non-universal amplitude ap- 
pearing in the description of bulk critical phenomena 42 -—, and 
P(d = 4) = 1/2 is a standard critical exponent. Within MFT 
T = t(^)^ 2 and u = 6a 2 (<i; l ~)~ 2 . In a finite-size system the 
bulk Hamiltonian J4?[<l>] is supplemented by appropriate sur- 
face and curvature (edge) contributions 42 - 43 . In the strong ad- 
sorption limit 49,50 , these contributions generate boundary con- 
ditions for the order parameter such that (j) 
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binary liquid mixtures these fixed-point (±) BC are the exper- 
imentally relevant ones. (Note that a weak adsorption prefer- 
ence might lead to a crossover between various kinds of effec- 
tive boundary conditions for the order parameter 7 ' 14 - 15 .) 

We have minimized numerically J4?[<j>] using a 3d finite 
element method in order to obtain the (spatially inhomoge- 
neous) profile m(r) for the geometries under consideration 
[see Fig.Q"). The normal and the lateral critical Casimir forces 
are calculated directly from these mean-field order parameter 
profiles using the stress tensor 30,35 . This allows one to infer 
the universal scaling functions of the critical Casimir forces at 
the upper critical dimension d = 4 up to an overall prefactor 
oc m _1 and up to logarithmic corrections. The corresponding 
critical Casimir potential is obtained by the appropriate inte- 
gration of the normal or of the lateral critical Casimir forces. 

In the case of planar walls the MFT scaling functions for 
the critical Casimir force can be determined analytically 35 
and one finds [see Eq. ©J for the case of symmetry break- 
ing boundary conditions the following critical Casimir ampli- 
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where K is the complete elliptic integral of the first kind, and 
A (+ % = -4A (+ +) [see Ref. [Hand Eq. (27) and Ref. [49] in 
Ref i El. 

In d — 4 (corresponding to MFT) the three-dimensional 
sphere is a hypercylinder and the physical properties are in- 
variant along the fourth dimension. Accordingly, the MFT 
results for the force and the potential given below are those 
per length along this additional direction. 



III. HOMOGENEOUS SUBSTRATE 

We first consider a three-dimensional sphere of radius R 
with (b) BC facing a chemically homogeneous substrate with 
(a) BC at a surface-to-surface distance D as shown in Fig.Q] 
denoting this combination by (a,b). The critical Casimir 
force F/ ab \(D,R,T) normal to the substrate surface and the 
corresponding critical Casimir potential <!>r a j,\(D,R,T) — 
JpdzFr ab \(z,R,T) take the scaling form a 6 ' 7 i 27 i 31 



R 



F m (D,R,T)=k B T^ s=T K ia j b) (e,A) (5) 



and 
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® {atb) (D,R,T)=k B T^Z2#(a, b )(®A), 



(6) 



where A = D/R and © = sign(f ) Dj%± (for t ^ 0) are the scal- 
ing variables corresponding to the distance D in units of the 
radius R of the colloid and of the correlation length £,±, respec- 
tively. The case d — 4 corresponds to the MFT solution up to 
logarithmic corrections, which we shall neglect here. Equa- 
tions (0 and © describe a force and an energy, respectively, 
per D d ~ 3 , which for d = 4 corresponds to considering F< a b \ 
and ®( a j)) P er length L4 of the extra translationally invariant 
direction of the hypercylinder. 



A. Derjaguin approximation 

The Derjaguin approximation (DA) is based on the idea of 
decomposing the surface of the spherical colloid into infinitely 
thin circular rings of radius p and area dS(p) = 2npdp which 
are parallel to the opposing substrate surface 6 ' 7 ' 27 ' 31 ' 51 . (Here 
we do not multiply 2npdp by the linear extension L4 of the 
hypercylinder along its axis in the fourth dimension, because 
the critical Casimir force is eventually expressed in units of 
La,, which therefore drops out from the final expressions.) The 
distance L of a ring with radius p from the substrate is given 
by 



L(p)=D+R il-Jl-p2/R2 



(7) 



Assuming additivity of the forces and neglecting edge effects, 
the normal critical Casimir forces dF(p) acting on these rings 
can be expressed in terms of the force acting on parallel plates 
[Eq. ©]: 



df(p) 
k R T 



dS 



[L(P)]' 



r *( 0>fc )(sign(f).L(p)/£fc). 



(8) 



Finally, in order to calculate the total force F{ a j>) acting on 
the colloid, one sums up the contributions of the rings, which 
yields 

F lah) (D,R,T) [R , 

{a ' h)y kBT L ^27tj o dpp[L(p)nV)(sign(f)L(p)/^ ± ). 

(9) 
(For d = 3, Fr a p, is the force on a sphere whereas in d = 4 it 
is the force on a hypercylinder per length of its axis.) 
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FIG. 2. (a) Scaling functions K/^ \ for the normal critical Casimir 
force [Eq. lO] acting on a three-dimensional sphere with (b) = (— ) 
BC close to a homogeneous substrate with (a) = (=p) BC IFig.fTl. 
The suitably normalized scaling functions K/^\ are shown as a 
function of the scaling variable = sign(f)Z)/£± fort ^ 0, where t is 

the reduced deviation from the critical temperature and Ki )(0,0) 

is the value of the critical Casimir force scaling function within the 
DA at T = T c for (— ,— ) BC. The solid lines correspond to the Der- 
jaguin approximation (DA, A = D/R — > 0) within mean-field theory 
(MFT, d = 4) whereas the dotted lines correspond to the DA obtained 
by using Monte Carlo (MC) results for films in d = 3 the systematic 
uncertainties of which are not indicated 52 . The normalization im- 
plies that at = both the solid and dotted lines pass through — 1 
for (— , — ) BC whereas the solid line passes through 4 for (+, — ) BC. 
The symbols correspond to the full numerical MFT results obtained 
for A = 1/3 and A = 1, the size of which indicates the estimated 
numerical error. (For (+, — ) BC and KOwe have not been able 
to calculate the corresponding scaling functions with adequate pre- 
cision due to severe numerical difficulties in obtaining the full three- 
dimensional order parameter profile in the presence of two "compet- 
ing" bulk values.) Since within the DA the dependence of Ki^- \ on 
A drops out, the difference between the symbols H and and the 
solid lines measures the accuracy of the DA in d = 4. (b) Differ- 
ence Ar> = iJ/ + _•) — ■&/ \ of the scaling functions for the Casimir 

potentials [Eq. l[6]l] for (+,— ) and (— ,— ) BC, suitably normalized 

by ■&/ )(0,0). The solid line corresponds to the DA within MFT 

and the symbols correspond to the full MFT results for A = 1/3 and 
A = 1 ; the dotted line is the DA for d = 3. Due to the normalization 
the solid line reaches 5 for = 0. 



One expects the DA to describe the actual behavior accu- 
rately if the colloid is very close to the substrate, i.e., for 
A = D/R —} 0. In this limit, Eq. (O can be approximated by 
L(p) = Da where a = 1 + p 2 /(2RD), so that one finds for 
the scaling function of the forc< 



,7,27 



K {aj7) (®,A^0)=27:J daa- f V,,)(a@), (10) 

and, accordingly, for the scaling function of the potential 67 - 31 



(11) 

At the bulk critical point, using Eq. (0, one finds 
the well known values Ki ab \(0,Q) — 2nA( ab \/{d—\) and 
d( a M(0,0) = 2itA( a>b -)/[(d-2)(d- 1)]. We note that the DA 
implies that the dependence of Fr a b \ and <£>t a h \ on the size 
R of the sphere reduces to the proportionality <* R indicated 
explicitly in Eqs. © and ©. 

Before proceeding further one first has to assess the accu- 
racy of the DA, which will carried out below within MFT 
(d — 4). We expect the range of validity of the DA to be 
similar for d = 3, so that within that range one can use the 
DA based on scaling functions for the film geometry obtained 
from Monte Carlo simulations 52 in order to calculate the crit- 
ical Casimir force acting on a colloid in d = 3. 



B. Scaling functions for the normal critical Casimir force and 
the potential 

The expressions obtained above within the DA hold for 
general boundary conditions (a) and (b) and are valid beyond 
the cases we consider in the following, i.e., a G { + , — } and 
b = — . Figure [2ja) shows the full numerical MFT (d = 4) 
results for the scaling functions ^(±,-) with A = 1 and 4 com- 
pared with the corresponding DA results based on the suit- 
able numerical integration [Eq. (TTUb l of the analytic (MFT) 
expression for h±--M. Moreover, in Fig. [2] the correspond- 
ing DA results for d = 3 are shown; they are obtained from 
the film scaling functions determined by MC simulations 52 
and by using the corresponding ratio of the correlation lengths 
above and below T,~. In Fig. |2jb) we report the difference 
A#(@,A) := i? (+ _)(@,A) - #j_ _)(€>, A) computed for the 
various cases reported in Fig. |2ja), which will be useful for 
describing the case of a chemically patterned substrate. The 
scaling functions in d = 4 are reasonably well reproduced by 
the DA for A < 0.4 and we expect this to hold for d = 3 as 
well. The fact that for increasing values of A the magnitude 
of the actual scaling functions becomes larger compared with 
those within the DA (corresponding to A — > 0) is in agreement 
with earlier results obtained for a ^/-dimensional hypersphere 
(see, e.g.,Ref.H3). 

It has been shown that the scaling functions obtained within 
the DA for d = 3 agree very well - within the experimental 
accuracy - with the ones obtained from direct measurements 
of the critical Casimir potential 6 -^ corresponding to A < 0.35 
(see also Ref. [48] in Ref. HI). 



IV. CHEMICAL STEP (s) 

The basic building block of a chemically patterned sub- 
strate of the type we consider here, i.e., with translational in- 
variance in all directions but one (x), is a chemical step (s) 
realized by a substrate with (a^) BC for x ^ at its surface. 
In this section we analyze the critical Casimir force if such 
a substrate is approached by a colloid with (b) BC with its 
center located at the lateral position x = X (see Fig. Q] and 
Ref. 32 for experimental realizations). We denote this config- 
uration by (a<\a>,b). The normal critical Casimir force F s is 
described by the scaling form 31 



R 



F t (X,D,R,T)=k B T^ s=T K t (S,&,A), 



(12) 



where S = X/\/RD is the scaling variable corresponding to 
the lateral position of the colloid. It is useful to write the scal- 
ing function K s as 



K S (S,&,A) 



K ( a< ,b) +K( a> ,b) 



K( a< ,b) —K( a> ,, 



-V(«k|«>.»)(S,®,A), (13) 



where the scaling functions of the laterally homogeneous sub- 
strates K( a> ,b) depend on © and A only [Eq. (|5j], and the scal- 
ing function W(a < \a > ,b) varies from +1 at S — > —°° to —1 at 
E — > +°o, such that the laterally homogeneous cases are recov- 
ered far from the step. Accordingly, the corresponding critical 
Casimir potential <$>s(X,D,R,T) = J£dzF s (X,z,R,T) can be 
cast in the form 3 -! 

<t> s (X,D,R,T)=k B T^A(Z,®,A), (14) 

and 



ft(S,0,A) 



fl( fl< ,6) + fl( a> ,6) 

2 
. %<,b) ~ #(«>,*) 



«(«<| a> ,i)(2,©,A), (15) 



where #( fl> M depend on © and A only [Eq. ©], and 
ffl( a< i a>] fe)(3 = ±°°,0,A) = =pl. Note that the scaling func- 
tions V / (a<|a>,6) an d (0 (a < \a > ,b) are independent of the common 
prefactortx uT x [see Sec. lIIBII . which is left undetermined by 
the analytical and numerical mean-field calculation of K s and 
A- 



A. Derjaguin approximation 

If the sphere is close to the substrate, i.e., A — >• 0, the DA 
can be applied, and one finds for the scaling function of the 
critical Casimir force [see Appendix lAl 

V( a< \ a> ,b)( z ^ 0,©,A^ 0) = Tl 

4/ 1 ~ S 2 /2 daa- rf arccos(|S|(2a-2)- 1 / 2 )M(a0) 



±- 



K {a<M) (©,A -> 0) -K {a>ib) (&,A -> 0) 



(16) 



where Ak(&) — k^ cl<b \(&) — fc( a> M(0) is the difference be- 
tween the scaling functions for the critical Casimir forces 
acting on two planar walls with (a < ,b) and with (a>,b) 
boundary conditions, respectively. We note that according 
to Eqs. ([Tol l and (TTOb within the DA W a< \ a> ,6) can be de- 
termined from the knowledge of the film scaling functions 
k(a,b)(®) [Eq. CJ] only. Due to the assumption of additiv- 
ity which underlies the DA, (i) yf( a< \ a> .b) vanishes at E = 
for all © and it is an antisymmetric function of E and (ii) 
W(a < \a > b) = W(a > \a < ,b)' within the DA both of these proper- 
ties are valid irrespective of the type of boundary conditions 
on both sides of the chemical step. (However, the actual scal- 
ing function W(a < \a > ,b) as > e -g-> obtained from full numerical 
MFT calculations may violate this symmetry because the ac- 
tual critical Casimir forces are non-additive.) At the bulk crit- 
ical point one has © = so that [see Appendix I A II . 

V / ( a< |« > ,fo)(S,© = 0,A^0) = 
E 2</-7 ( is (3 _ d) + (3 _ 2d )Z 2 -Z 4 ) (2 + Z 2 y {d ^- ] (17) 

independent of k< a> b y Similarly, within the DA one finds for 
the scaling function CO of the critical Casimir potential [see 
Appendix[A]and Ref. 13111 



fl>(o<|a>,6)(S ^ O,0,A^ 0) = Tl 



±- 



"4 f°° a 5arcc08 ( 5 ' ) — \A— 1 



Ait(©[l+EV 2 ] 



i? (fl<)fe) (0,A^O)-t? (a>>6) (©,A^O) 



(18) 



This yields C0i u< \ a> m(E = 0,©,A — > 0) =0, as expected from 
the underlying assumption of additivity; within full MFT this 
only holds in the limit A — ► 0. At the critical point we find [see 
Appendix I A II 



CO, 



(«< 



fl>ifc) (E,0 = O,A->O) =Z(l-d-Z 2 )(Z 2 + 2 



-3/2 

(19) 



For symmetry breaking (T, — ) BC and S> 1 the critical 
Casimir force ft^ _\ (D, T) acting on two planar walls at a dis- 
tance D decays °c exp(— ©) [Eqs. (0 and Q)L which within 
the DA leads to the same d-independent result for the scaling 
functions W+\ ) and &u.\ ) [see Appendix lA2l : 

Vr (+ |__)(S,0»1,A->O) = 

0) (+ |_ _)(S,0 >1, A ->0) = 



-erf V0/2S 



(20) 



where erf is the error function. 

Figure|3ja) compares the scaling function CO^ a< \ a> ,b) for the 
critical Casimir potential of a sphere with (— ) BC in front of 
a (+|— ) step, as obtained within the DA for d = 4 [Eq. (fT8bl. 
with the one obtained numerically within full MFT for A = 
1/3. For A < 1/3 the DA captures the scaling function very 
well, in particular for © > 3 3 ' . The scaling function C0r a< \ a> b \ 
obtained within the DA (d = 3) on the basis of the Monte 
Carlo data of Ref. |2j, which is also shown in Fig. |3ja), has 
been used successfully in order to interpret the experimental 
data of Ref. |32j, for which the analysis in terms of separate, 
independent, and consecutive chemical steps turned out to be 
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FIG. 3. (a) Scaling function ft)/ + i \ [Eq. J15b l for the critical 

Casimir potential of a spherical colloid with (— ) BC across a chem- 
ical step (+|— ) as a function of S = X/\/RD for various (positive) 

values of = D/E,^2L. Within the DA (Ot + i \ is an antisymmetric 

function of E [Eq. dl8H whereas within full MFT this antisymme- 
try is slightly violated, in particular for small 0. (b) Corresponding 
scaling function K$ [Eq. J21M of the lateral critical Casimir force, 

normalized by the amplitude Ki )(0,0) = 27rA( \/{d — 1) of 

the normal critical Casimir force atf = T c acting on a colloid with 
(— ) BC close to a homogeneous substrate with (— ) BC within the 
DA [Sec. HITAl . For both (a) and (b) the full numerical MFT re- 
sults obtained for A = 1/3 are shown as symbols (the symbol size 
represents the estimated numerical error) whereas the lines show 
the corresponding results obtained within the DA (i.e., A — > 0); the 
dotted lines refer to d = 3 and are obtained by using Monte Carlo 
simulation data 5 and the solid lines refer to d = 4. The lines for 
= are obtained by using Eq. d 19b and Eq. ( |23t , respectively; for 
= 3.2,4.7,8.1 the DA lines de facto coincide with the asymptotic 
results obtained for symmetry breaking BC and 0^1 [Eq. l |20| l and 
Eq. d24t . respectively] and thus are indeed independent of d. The DA 
(d = 4) provides a good approximation for the full numerical MFT 

data, in particular for > 3. K^ > implies that the colloid moves to 
the right where it enjoys an attractive potential versus a repulsive one 

for H < 0. Within the DA IQ is a symmetric function of S [Eqs. j!8t 
and d22H whereas within full MFT this symmetry is slightly violated, 
in particular for small 0. 



accurate. Moreover, the critical Casimir forces turned out to 
be a sensitive probe of the chemical pattern and its geometric 
design 31 . 



B. Lateral critical Casimir force 

The lateral critical Casimir force is given by F s " = — dx&s 
and can be cast in the scaling form 



R ID 



1/2 

F?(X,D,R,T)=k B T-^ rT l-) ^'(3,0, A), (21) 



where K's is a universal scaling function. F s " and A^ 1 vanish 
far from the chemical step, i.e., for |S| — > °°. In Eq. d2Tb the 
prefactors in terms of R and D and their exponents are chosen 
such that Kg is regular and non-vanishing for A — > 0. We note 
that the same holds for the normal critical Casimir forces and 
the corresponding potentials [see Eqs. (0, ©, (fT2l) . (fT4l >. and 
the considerations following below]. 

Within the DA Kjj can be calculated from Eqs. ( TT~5b and 

tf|(E,0,A->O) = 

- \ [#(«<,*) (®, A -»• 0) - $ {a>yb) (0, A -> 0)] x 

^ffl (a< | a>i6) (S,0,A^O). (22) 

At bulk criticality © = one finds with Eq. dT~9T > [see Eq. ([TBT ll 



£j(S,0 = O,A->-O) = 7rM(O)(2 + E 2 ) ( '' ^ 



(23) 



For (=F, -) BC and > 1 Eqs. (O, (fB), and © lead to 



#J(S,©>1, A -)■()) = 

[t? (+] _ ) (0,A)-t? ( _ ) _ ) (0,A) 



exp 



&Z 2 



(24) 



for /?of/i d = 3 and d = 4. [The pref actor A#(0,A) = 
%-)(©, A) - t> ( _ _)(0,A) in Eq. © is shown in Fig.Hb).] 
Figure Ob) shows the comparison between the normalized 
lateral critical Casimir force obtained within the DA (solid 
lines) and the full MFT data obtained for A = 1/3 (symbols). 

We infer that not only the shape of K$ as a function of S but 
also its amplitude is described well by the DA [Eqs. d23l and 
(1241 1 for A < 1/3, and in particular for > 3. We expect this 
feature to hold in d = 3, too, as well as for the normal critical 
Casimir force and the critical Casimir potential. The lateral 
critical Casimir forces for d = 3 obtained within the DA on the 
basis of Monte Carlo simulation data for the film geometry 52 
are shown in Fig. [3jb) as dashed lines. Compared with the 
previous curves, these ones have similar shapes but their over- 
all amplitudes in units of the normal critical Casimir force at 
= are significantly different for = and = 3.2. This 
difference reflects the analogous one observed in the normal- 
ized difference between the corresponding critical Casimir po- 
tentials for (+, — ) and (— ,— ) BC, reported in Fig.[2|b). 



V. SINGLE CHEMICAL LANE (() 

In this section we consider the case of a colloid with (b) BC 
close to a substrate with a single chemical lane (£) with (a() 
BC and width 2L in the lateral x direction and which is invari- 
ant along the other lateral direction(s). The remaining parts of 
the substrate are two semi-infinite planes at \x\ > L with (a) 
BC [see Fig.|T}. The lateral coordinated of the center of mass 
of the sphere along the x direction is chosen to vanish in the 
center of the chemical lane. One expects that for "broad" lanes 
a description in terms of two subsequent chemical steps is ap- 
propriate [Sec.lIVIand Ref. |3l||, whereas for "narrow" lanes 
the effects of the two subsequent chemical steps interfere. We 
find that in addition to the variables characterizing the chem- 
ical step [Eq. ([T2ll. a further scaling variable A = L/^fRD 
emerges naturally, which corresponds to the width of the lane. 
Accordingly, the normal critical Casimir force Fi acting on the 
colloid can be cast in the form 

F e (L,X,D,R,T)=k B T-^ T K ( (A,Z,©,A), (25) 

where Kf is the corresponding universal scaling function. The 
critical Casimir potential scales as 



R 



® e (L,X,D,R,T)=k B T^ [ZI #e(A,Z,®,A), 



(26) 



with fy as the universal scaling function for the potential of 
a sphere close to a single chemical lane. Analogously to 
Eqs. ( fT3] l and ( fT5l ) we define yf£ an d ®e according to 



^(A,S,0,A) = 



K(a.b) +K(a e ,b) 



K {a,b) - K (a t ,b) 



%(A,S,0,A), (27) 



and 



a*(A,s,e,A) = 



%,b) + % e ,b) 



®(a,b) - &(a e ,b) 



a*(A,E,®,A), (28) 



so that far from the lane Ye{A, |S| ^> A,0,A) = 
0%(A, |S| 3> A,0,A) = 1. On the other hand, only 
for a "broad" lane the scaling functions at the cen- 
ter of the chemical lane approach their limiting value 
Yi(A-K»,S=O,0,A) = -1 = ffl f (A^oc,S = O,0,A), 
corresponding to the homogeneous case with (d(,b) BC. 



A. Derjaguin approximation 

Using the underlying assumption of additivity of the forces, 
within the DA (A — > 0) we find for the scaling functions of the 
critical Casimir force and of the critical Casimir potential [see 
AppendixlBl 

VG(A,S,®,A->0) = 

1 + V(a f \a.b) (S + A, 0, A -> 0) - W(ae\a,b) (S - A, 0, A -»■ 0) 

(29) 



and 



l+®(a c \aJ,)(Z+A,©,A^O)-co iaf \ a . h) (Z-A,®,A^O), 

(30) 

respectively. Thus, within the DA, from the knowledge of the 
scaling functions W(a t \a,b) [Eq. ( Q ] and (0t ai \ a ,b) [Eq- dHl >] for 
the chemical step with the appropriate BC, one can directly 
calculate the corresponding scaling functions for the chemical 
lane configuration. Accordingly, in the limit A — >• and for 
symmetry breaking BC, y^ and (0( can be analytically calcu- 
lated on the basis of Eqs. J29l ) and ( |30l by taking advantage of 
Eqs. {FT), ([19), and (EDJ. 



B. Scaling function for the critical Casimir potential 

In Fig. Ua) we show the scaling function C0( for the crit- 
ical Casimir potential obtained within the DA for d = 3 and 
d = 4 (MFT) at the bulk critical point T = T c [Eqs. (f30) and 
( TT9l l1 for various values of A = Lj \JKD as a function of the 
lateral coordinate of the colloid. One can infer from Fig. [4] 
that, at bulk criticality, the critical Casimir potential varies 
less pronounced in d = 3 than in d = 4. As expected, for 
small values of A (i.e., "narrow" chemical lanes), the poten- 
tial does not reach the limiting homogeneous value — 1 in the 
center of the chemical lane. On the other hand for large val- 
ues of the scaling variable A (i.e., "broad" chemical lanes), 
CO( does attain the value — 1 in the center of the chemical lane 
and the critical Casimir potential flattens. In this case the po- 
tential is adequately described by two independent chemical 
steps. However, the criterion for being a sufficiently "broad" 
lane depends sensitively on © and d. Indeed, from Eqs. ( 1301 
and $1% we find that at criticality (© = 0) the critical Casimir 
potential at the center of the chemical lane (E = 0) reaches 
the limiting value corresponding to the colloid facing a homo- 
geneous substrate by up to 1% for A>3.3inJ = 4 and for 
A > 10 in of = 3. We note that the curves in Fig. HJa) as well 
as these bounds are independent of the actual boundary con- 
ditions because for all kinds of BC the scaling function of the 
normal critical Casimir force is constant at the critical point 
[see Eq. (fJJ]. 

Below we shall discuss some properties which are specific 
for BC with a,a£,b £ {+,—}, which exhibit the feature that 
the normal critical Casimir force /fe -) acting on two pla- 
nar walls decays purely exponentially [see the text preceding 
Eq. (01] as a function of their distance expressed in units of the 
bulk correlation length [see Eqs. (fll and (J3j] . In Fig.[4fb) the 
scaling functions G)g in d — 3 and d = 4 obtained from Monte 
Carlo simulation data 52 and analytic MFT results 35 , respec- 
tively, within the DA [see Eqs. (1301 and ( 1201 )1 are shown for 
the same values of A as in Fig. H|a) but off criticality. For 
© = 7.7 the curves for d — 3 and d = 4 are indistinguishable 
from each other and from their common asymptotic expres- 
sion given in Eq. ( 1201 [see also Ref. 31] . For © 3> 1, the crit- 
ical Casimir potential attains its limiting homogeneous value 
in the center of the lane for values of A which are smaller 
than the ones for = due to the shorter range of the forces. 
That is, for both d = 3 and d = 4 the single chemical lane is al- 
most equally well approximated by two independent chemical 
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FIG. 4. Scaling function (Of> [Eq. d28H describing the lateral variation 
of the critical Casimir potential of a colloid across a single chemical 
lane of width 2L as a function of the lateral position X of the col- 
loid in units of the half width of the lane [see Fig.Q] S = X/\/RD, 
A = Lj \JRD, = £>/£+]. Here, cog has been obtained within the 
DA (A — > 0) in d = 3 and 4 [Eq. J30H. In (a) the curves correspond 
to = [Eq. dl9H, whereas in (b) they correspond to = 7.7 and 
a,a e ,b 6 {+,-} BC [Fig.Q]. For > 1 [(b)] the corresponding 
scaling functions obtained from Monte Carlo simulation data 52 in 
d — 3 and from analytic MFT results in d — 4 de facto coincide 
and their asymptotic expressions are given by Eqs. l !20t and d301 >. 
(£)( = 1 corresponds to the laterally homogeneous critical Casimir 
potential for (a, b) BC outside the chemical lane, whereas C0( = — 1 
corresponds to the value of the critical Casimir potential for the ho- 
mogeneous case with (ag,b) BC as within the chemical lane. For 
large values of A the critical Casimir potential is the same as for two 
independent chemical steps, and a>g reaches its limiting value — 1 in 
the center of the lane at E = [see the main text]. In (b), for 0>1, 
a>( attains — 1 in the center of the chemical lane already for smaller 
values of A due to the exponential decay of the critical Casimir force. 
We note that the DA results for = (i.e., at the critical point) are 
independent of the actual boundary conditions which, accordingly, 
were not specified in (a). 



steps for A > 1 .5 at = 3.3 (data not shown) and for A > 1 .0 
at = 7.7 [Fig. Hb)]. 

In Fig. [5] we compare the MFT C0( obtained within the DA 
(A -» 0) at © = [Eqs. (f30]l and (Q30l] with the scaling func- 
tion obtained from the full numerical MFT calculations for 
A = 1/3. We find a rather good agreement even for small val- 
ues of A (i.e., "narrow" chemical lanes). This shows that for 
the geometry of a colloid close to a single chemical lane, non- 
linearities, which are actually present in the critical Casimir 
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FIG. 5. Test of the performance of the DA for the scaling function Cty 
[Eq. J28H of the critical Casimir potential for a sphere with (— ) BC 
close to a single chemical lane with (— ) BC embedded in a substrate 
with (+) BC. The MFT m t is evaluated at bulk criticality = in 
d = 4 both on the basis of the DA (lines, A — > 0) and of the full numer- 
ical MFT (symbols, A = 1/3). There is good agreement between the 
DA and the full MFT results, even for small values of A = L/\/RD. 
Nonlinear effects, which are inherently present in the theory, do not 
strongly affect the potential. For A — > the assumption of additivity 
of the critical Casimir forces underlying the DA is reliable even for 
small A. 



effect and potentially invalidate the assumption of additivity 
underlying the DA, do not affect the resulting potential for 
small values of A. We expect this property to hold beyond 
MFT in d = 3 as well, in particular off criticality, i.e., for 
© ^ 0. 



VI. PERIODIC CHEMICAL PATTERNS (p) 

In this section we consider a pattern of chemical stripes 
which are alternating periodically along the x direction. The 
pattern consists of stripes of width L\ with (a\) BC joined 
with stripes of width Li with (02) BC, such that the periodic- 
ity is given by P = L\ +L2- Thus, the geometry of the sub- 
strate pattern is characterized by the two variables L\ and P 
[see Fig.Q"). The coordinate system is chosen such that the 
lateral coordinate X of the center of the sphere is zero at the 
center of a (ai) stripe. The normal critical Casimir force F p 
acting on the colloidal particle and its corresponding potential 
<I>p take on the following scaling forms: 

F p (L h P,X,D,R,T) =k B T^ ZI K p {X,n,S,®,A) (31) 



and 



%(L U P,X,D,R,T) =& B r^2# p (A,n,S,0,A), (32) 

where II = P/^fRD is the scaling variable characterizing the 
periodicity of the pattern and A = L\jP is the scaling vari- 
able chosen to correspond to the relative width of the stripe 
with («i) BC. K v and # p are universal scaling functions for 
the normal critical Casimir force and the critical Casimir po- 
tential, respectively. For A — 1 or the force and the poten- 
tial correspond to the homogeneous cases with {a\,b) BC or 



(ci2,b) BC, respectively [see Sec. HTTll . As before it is useful to 
define scaling functions \f/ v and C0p which vary for A € [0, 1] 
within the range [—1,1] and describe the lateral behavior of 
the critical Casimir effect: 



Xp(A,n,S,0,A) 



K (a 2 ,b)+ K { ai ,b) 



K (a 2 ,b) ~ K {a u b) 



VAp(A,n,S,©,A) (33) 



and 



# p (A,n,S,0,A) = 



0(a2,&) + &{ai]b) 



®{a 2 ,b)-®{a h b) 



fflp(A,n,S,0,A). (34) 



A. Derjaguin approximation 

Taking advantage of the assumption of additivity of the 
forces underlying the DA, one finds for the scaling function 
of the normal critical Casimir force in the limit A-)0 [see 
AppendixICl 

V/ p (A,n,S,0,A^O) = 

l + E {^( fll Kfe)(S+n(n + |),0,A^O) 

n=— <» 

-^( flI |a 2 ,i.)(S + n(n-|),0,A^O)}. (35) 

Thus, the knowledge of the scaling function ty( ai \ a2 fi ) for a 
single chemical step with the appropriate BC [Sec. lIVI is suf- 
ficient to calculate directly the corresponding scaling func- 
tion of the critical Casimir force acting on a colloid close 
to a periodic pattern of chemical stripes. As expected, from 
Eq. ( f35T > one recovers the values y/ p (A =O,n,E,0,A) = 1 
and y/n(A = l,n,E,0,A) = — 1, i.e., the cases of a colloid 
with (b) BC facing a homogeneous substrate with (02) BC 
and (a\) BC, respectively [see AppendixICl. 

In the limit II — > 0, i.e., for a pattern with a very fine struc- 
ture compared to the size of the colloid, the sum in Eq. d35l ) 
turns into an integral [see AppendixICl and, as expected, y/ p 
becomes independent of S, i.e., of the lateral position of the 
colloid: 



V/ p (A,n^O,E,0,A^O) = l-2A. 



(36) 



Accordingly, in the limit II — >• the force acting on the colloid 
- within the DA - is the average of the ones corresponding to 
the two boundary conditions weighted by the corresponding 
relative stripe width [see Eqs. (l36l l and (l33ll: 

/: p (A,n^o,s,©,A^o) = 

Li K {aub) (&,A^0) + — ^-# (a2i6) (0,A^O). (37) 



L1+L2 



L1+L2 



For the scaling function of the critical Casimir potential the 
results are completely analogous to Eqs. d35ll-(l37]i [see Ap- 
pendixlcl. 
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FIG. 6. MFT (rf = 4) scaling function i|/ p [Eq. OS] of the normal 
critical Casimir force acting on a colloidal sphere with (b) = (— ) 
BC which is close to a periodically patterned substrate [Fig.fT) with 
(ai) = (— ) BC on one kind of stripes [shaded areas] and (02) = (+) 
BC on the other kind of stripes. Due to this choice of the BC the 
colloid is attracted by the shaded stripes and repelled by the others. 
\ffp is shown as a function of the lateral position of the colloid X/P 
with P = L\ + Li and at the bulk critical point = 0. The geom- 
etry of the pattern is characterized by n = P/\/RD and X = L\ /P, 
for which we have chosen the values (a) X = 0.5 and (b) X = 0.2. 
The lines are the results for \j/ p as obtained within the DA for d = 4 
[Eqs. d35t and dl7H . whereas the symbols represent the full numer- 
ical data obtained within MFT for A = 1/3 for various values of IT. 
For patterns which are finely structured on the scale of the colloid 
size, i.e., IT < 2, the actual results deviate from the approximate ones 
obtained within the DA due to the strong influence (in this context) 
of the inherent nonlinear effects. 



B. Scaling function for the normal critical Casimir force 

Figure[6]shows the scaling function y/ p [Eq. (l33l ll as a func- 
tion of S/II = X/P, describing the lateral variation of the 
normal critical Casimir force at = as obtained within the 
DA for d = 4 [Eq. ( f35T > with Eq. ([PTT i; solid lines] compared 
with the one obtained from the full numerical MFT calculation 
[A = 1/3; symbols] for symmetry breaking boundary condi- 
tions (01) = (-), (a 2 ) = (+), and (b) = (-) [Fig.Q]. From 
this comparison for A = 0.5 [Fig.[6]a)] and A = 0.2 [Fig.^b)] 
and for various values of IT one can infer that for A — > and 
IT > 1, i.e., Lx+L 2 ^> VRD the DA describes well the ac- 
tual behavior of the scaling function, even if the force scal- 
ing function does not attain its limiting homogeneous values 
y/p = ±1 in the center of the stripes. However, for IT < 2 
(in d = 4 at T = T c ) the DA does not quantitatively describe 
the actual behavior and the scaling function y/ p obtained from 
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FIG. 7. (a) The same as in Fig. [6] but for X = 0.8. Also in this 
case, the DA turns out to be accurate for IT > 2 while it fails to de- 
scribe quantitatively the full numerical data for smaller values of IT. 
(b) Comparison between the scaling functions i// p in d = 3 (dotted 
lines) and d = 4 (solid lines), at T = T c , for X = 0.8, and within the 
DA. At the critical point the expression for this scaling function l//p 
is known analytically [see Eqs. d35b and JT7H. and the corresponding 
plot presented here shows that the lateral variation of the normal crit- 
ical Casimir force is less pronounced in d = 3 than in d = 4. (We note 
that for IT — > we expect that also in d = 3 the DA fails to describe 
quantitatively the actual behavior; however, we nonetheless present 
the curve for IT = 0.57 in order to show that the critical Casimir force 
obtained within the DA practically does not change laterally for such 
small values of IT.) 



the full numerical MFT calculations deviates from the one ob- 
tained within the DA. Within both the DA and the full numer- 
ical MFT calculation, for II — > the normal critical Casimir 
force loses its lateral dependence on S. But from the full nu- 
merical calculation we find that the corresponding constant 
value which is attained by y/ p differs from the one obtained 
within DA [Eq. (I36il1. This shows that for small periodicities 
P < \/RD nonlinearities inherent in the critical Casimir effect 
strongly affect the resulting scaling functions of the force and 
the potential, so that in this respect the assumption of additiv- 
ity of the force and thus the use of the DA are not justified. 

Figure [TJ a) shows the same comparison as Fig. [6] but for 
A = 0.8, which corresponds to an areal occupation of 80% of 
the substrate surface with (-) BC and 20% with (+) BC. Due 
to the fact that at the critical point w ai i a2b \ (S, © = 0, A — >• 0) 
is actually independent of the BC, i^ p (A = 0.8,IT,E,© = 
0, A — >• 0) in Fig. |7Ja) is, within the DA, complementary to 
the one for A = 0.2 in Fig. |6jb), i.e., it is obtained from the 
latter by a reflection with respect to y/ p = followed by a shift 
in S/II of 0.5. Instead, the full numerical data in Fig.|7|a) and 
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FIG. 8. Normalized scaling function A 1 of the critical Casimir 
force at criticality acting on a homogeneous planar wall with (+) 
BC opposite to a periodically patterned planar substrate with stripes 
of alternating (+) and (— ) BC as a function of v = 5-/5+, where 
5+ and 5- are the respective widths. The symbols correspond to the 
MFT (d = 4) data presented in Fig. 12 of Ref. 33 for various values 
of 5+/L (note that A ( |+ = A (++) /(d - 1) in Fig. 12 of Ref.H3). The 
dashed and dotted lines which join the data points are a guide to the 
eye. The solid line corresponds to the DA result given in Eq. J41t 
which assumes additivity of the forces and turns out to be indepen- 
dent of the ratio 5+ /L. One can immediately infer from the graph 
that here the assumption of additivity is not justified, which is the 
limiting configuration of the sphere-wall geometry for IT — > 0. 



Fig- 13b) show a different behavior as they clearly tend to as- 
sume the value — 1 corresponding to the homogeneous case 
with (— ,— ) BC. By contrast, for the case A = 0.2 shown in 
Fig. |6fb), the full numerical data do not reach as closely the 
value +1 corresponding to (+, — ) BC, although the substrate 
area is covered by 80% with (+) BC. This feature is addressed 
in more detail in Sec. IVIII Figure [T^b) compares the scaling 
function i// p of the normal critical Casimir force at T — T c and 
for A = 0.2 as obtained within the DA for d = 4 (solid lines) 
with the corresponding one for d = 3 (dotted lines). At T = T c , 
\jfp is determined by Eqs. fl35l > and ( fTTl l from which one can in- 
fer that the lateral variation of the normal Casimir force is less 
pronounced for d = 3 than for d = 4. This qualitative feature 
holds for all values of A (not shown). However, off critical- 
ity, © > 1, [according to Eqs. ([35j and <j20j>] the DA scal- 
ing functions both for d = 3 as obtained from MC simulation 
data and for d = 4 as obtained from MFT de facto coincide 
(not shown), similarly to the case of a single chemical lane in 
Fig. Kb). 

Although one would expect the DA to be valid for large 
radii R, the lateral variation of the boundary conditions at the 
surface of the patterned substrate on a scale P < \/RD - corre- 
sponding to the limit TI — > - renders the DA less accurate, as 
it clearly emerges from the numerical data presented in Figs. [6] 
and|7] The fact that a large colloid radius R does not guarantee 
the validity of the DA can be understood by noting that such 
a discrepancy between the full numerical calculation and the 
result of the DA approximation already emerges in the film ge- 
ometry (formally corresponding to the limit R — > oo), i.e., for a 



chemically/patterned wall opposite to a laterally /zomogeneous 
flat wall. This "ph" configuration has been studied in Ref.uJ 
within MFT for laterally alternating chemical stripes of width 
L\ = S+ and L^ = S- with (+) and (— ) BC, respectively, op- 
posite to a homogeneous substrate with (+) BC a distance L 
apart [see Fig.Q]and the inset of Fig. [8). Indeed, by using the 
assumption of additivity of the critical Casimir forces under- 
lying the DA and neglecting edge effects, the normal critical 
Casimir force fjL.JS+,S-,L,T) per unit area acting on the 
walls is predicted to be given by 



/(DA)(^+'^-'^' r ) 



-s_ 



/(+,+) (L, T)- 



S- 



S++S- 



/(+,_) (L,T), (38) 



where /(+,±) refer to homogeneous parallel walls, as in 
Eq. (HJ. At the bulk critical point the critical Casimir force 
is given in general by22 

f h (S + ,S-,L,T = T c ) = k B T c ^-A»" (v = ^-, ?± 

' ' (39) 

Using Eq. ( f38l > together with Eqs. (Q]) and (O one finds within 
the DA that 



(d-l)A' 



P h 

(DA) 



vA (+,-)+ A (H 
1+V 



(40) 



which renders the rhs of Eq. ( f40l > to be independent of the 
scaling variable S+/L. Within MFT as studied in Ref.[33|(a? = 
4), one has A( + \ = — 4A( + + ) > [see the end of Sec. Ill B II 
so that 



A ^L 



|A(+.+)| 4v-l 



(41) 



In Fig. [8] we show the comparison between the actual scaling 
function A ph (data points, obtained numerically as reported in 
Fig. 12 of Ref. HI) and A^ A) (Eq. (gUl, solid line) derived by 
assuming additivity of the forces and neglecting edge effects. 
Figure [8] clearly shows that the actual behavior of the critical 
Casimir force in the film geometry is not properly predicted 
within these assumptions. This is expected to be due to the 
presence of nonlinear effects and of edge effects in this con- 
text. This explains why in the limit TI ->• the DA (R > D) 
used here does not capture the behavior of the critical Casimir 
force acting on a colloid close to periodically patterned sub- 
strate. 

In Fig. [9] we show the behavior of scaling function K p 
[Eq. (13TT >1 of the normal critical Casimir force acting on the 
colloid in d = 4 with (b) = (-) BC as a function of© = D/| + 
(i.e., as a function of the normal distance of the colloid from 
the substrate in units of the bulk correlation length) and for 
various values of A and II, In Fig.|9]the scaling function K p is 
evaluated at X = [see Fig.[T| which corresponds to the most 
preferred lateral position of the colloid in which the normal 
force is least repulsive or most attractive [see Fig. [6]]. From 
Fig. [9] one can infer that the DA does not provide an accurate 
estimate of K p in the whole range of for TI = 0.57 [panel 
(b)], whereas it does so for TI = 2.3 [panel (a)]. Indeed, for 
IT = 0.57 the discrepancy between the DA and the numer- 
ical data is already significant for < 4 and 0.3 < A < 0.9, 
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n = 2.3 lines: DA (A->0) (a) 

: = symbols: MFT (A = ±) 
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■4 = 2/5 
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FIG. 9. Scaling function &" p [Eq. J31H of the normal critical Casimir 
force acting on a spherical colloid with (— ) BC located at X = 
(E = X/\/RD) close to a periodically chemically patterned substrate 
[see Fig. [T). K v is suitably normalized by the absolute value of the 
force scaling function Ki )(0,0) = 2:z/S.i\/{d — 1) for the ho- 
mogeneous (— , — ) case at criticality and within the DA [Sec. lIII A"l . 
The lateral position of the center of the colloid is fixed at the center 
of a stripe with (a\) = (— ) BC and width L\ = XP, which it is at- 
tracted to, in contrast to the second type of stripes with (02) = (+) 
BC and width Lq, = (1 — X)P, which it is repelled from. The scaling 
variable corresponding to the periodicity of the substrate pattern is 

(a) IT = P/\/RD = 2.7 and (b) IT = 0.57, whereas the relative area 
fraction of the (— ) stripes changes from X = L\/(L\ +L2) = to 
1 = 1 (top to bottom: fully repulsive to fully attractive). In (a) and 

(b) the lines represent the result for the MFT critical Casimir force 
within the DA [A — > 0, d = 4, see Eq. d35H, whereas the symbols 
represent the full numerical MFT data obtained for A = 1/3. The 
DA agrees reasonably well with the full data for Ft = 2.3 [(a)] and 
© > 1, but for Ft = 0.57 [(b)] it fails to describe the actual behav- 
ior within the ranges < 4 and 0.3 < X < 0.9 where the nonlinear 
effects strongly affect the resulting scaling function. In (c) K p is 
shown for TI = 0.57 and 2.3, as obtained for d = 3 within the DA on 
the basis of the Monte Carlo simulation data for the film geometry 52 . 
(We note, however, that we do not expect that the curves shown for 
IT = 0.57 are quantitatively reliable.) 



whereas for TI = 2.3 agreement is found for all values of X ex- 
cept for © < 1 [Fig.|9ja)]. This fact suggests that for relatively 
small periodicities TI < 2 non-additive and edge effects be- 
come important. On the other hand, for large values of © 3> 1 
the DA describes the behavior of K p rather well for all values 
of TI due to the exponential decay of the critical Casimir force 
for 0^1 [Eq. (0]. Figure[9fc) shows the scaling function K p 
for d = 3 within the DA as obtained from Monte Carlo simu- 
lation data for the film geometry 52 . The qualitative features of 
the behavior of K p in d — 3 and d — 4 are similar. 

From our analysis in d — 4 we conclude that the DA de- 
scribes quantitatively well the behavior of the actual critical 
Casimir force for TI > 2 for all values of ©. For smaller val- 
ues of TI, the DA is only quantitatively reliable for large values 
of © (at which the force decays exponentially). For example, 
for TI > 0.5 the DA result is quantitatively correct for > 4. 
We expect these properties to be carried over to d = 3. 



C. Critical Casimir levitation 

Rather remarkably, within a certain range of values of X, 
Kp changes sign as a function of © = D/4+ [Fig. |9). In this 
context it is convenient to introduce for later purposes another 
scaling variable ^ — Til©! 1 ' 2 = P j \JR^± which is indepen- 
dent of D and therefore does not vanish in the DA limit D^lR 
(i.e., A — > 0). Due to this change of sign of K p , there exists a 
certain value = ©oOP,A,S,A) at which the normal critical 
Casimir force F p acting on the colloid vanishes. This implies 
that in the absence of additional forces the colloid levitates at 
a height Do determined by ©0 and <jj + , which can be tuned 
by changing the temperature. Since for fixed geometrical pa- 
rameters R, X, and P the scaling variables 0, II, S, and A 
depend on D, one has to consider the behavior of F p as a func- 
tion of D near Do in order to assess whether the levitation 
is stable against perturbations of D or not. Stability requires 
doFp \d=d < (so that for D <Dq the colloid is repelled from 
the patterned substrate, whereas for forD > Do it is attracted). 
According to Eq. OTb one has 



R 



d D F p =k B Tj^x{-(d-l) 

-±Hdn-±Zd z + ®d @ + Ad A }K p {X,Tl,Z,®,A). (42) 

The laterally preferred position is always at X = Xq = 0, cor- 
responding to E = So = 0, so that within the DA (A — » 0) one 
has 



sign 



^=D Q J(=X ,DA 



[d D F p \ D _ 
gn({-£n£i + 0fc}tfp(A > II > E = O,0 > A->.O)| e= 



0„ 



(43) 



where we have used the implicit equation F p \o=d q = so that 
^pId=d ~ 0. (Equation (l43l assumes that d&K p does not di- 
verge oc A -1 for A — > 0.) In the following we only consider 
> and BC (at) = (-), (a 2 ) = (+), and (b) = (-). 

Within the DA we find that both <3nXp|© = © 0j s=So an d 
<50^p|©=@ o ,s=s o are negative, so that according to Eq. d43l the 
sign of doF p \D=D ,x=x ,DA can vary and depends on their val- 
ues as well as on ©o and II, However, at criticality (0 = 0) the 
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FIG. 10. Values of the scaling variable ©0 at which within the DA 
(A — > 0) the normal critical Casimir force A" p shown in Fig. [9] van- 
ishes as a function of Y for (a) d = 4 and (b) d = 3 on the basis of 
Monte Carlo simulation data 52 and for various values of X — L\/P. 
The solid lines correspond to values of ©o for which the levitation 
of the colloid at a height Dq above the substrate is stable against 
perturbations of D [doF ? \o=D < 0* see Eq. J43H , The shaded re- 
gion and the dashed lines indicate those values of ©0 for which 
<?d^pId=Z5 > ar, d tnus do not correspond to stable levitation. For 
X > Xq with Xo(d = 4) = 4/5 and Ao(rf = 3) ~ 0.88, ® ceases to ex- 
ist, i.e., K ? does not exhibit a zero. For X < X\ with X\ (d = 4) = 1/2 
and Xi{d = 3) ~ 0.545, 0o(*P \ *o(A)) diverges. (The values for 
'i'rj(A) are indicated by upward arrows.) For any X < Xq, ©0 ex- 
ists for ¥ < y ¥*(X). (From the analysis in Fig. [5] we expect the DA 
to be quantitatively reliable only for *J? > 2y / ©o for ©o < 4 and for 
* > 0.5 V©o for ©o > 4, which implies X < 0.7 in d = 3 and X < 0.6 
inrf = 4.) 



second term of the rhs of Eq. d43l vanishes. Thus, at the bulk 
critical point T = T c the derivative dpFp evaluated at D = Dq 
and X = Xq = is always positive so that one cannot achieve 
stable levitation. On the other hand, for © > it is always 
possible to find geometrical configurations for which the col- 
loid exhibits stable levitation, as described in the following. 



Figure[l0] shows the values of ©0 at which the normal criti- 
cal Casimir force acting on a colloid vanishes as a function 
of the new scaling variable \P introduced at the beginning 
of this subsection, for various A, for E = 0, and within the 
DA (A -> 0) for (a) d = 4 and (b) d = 3. The correspond- 
ing sign of doF p \ D=D [according to Eq. d43l l is also indi- 
cated: @o drawn as a solid line indicates doF v \ D=D < 0, 
i.e., stable levitation of the colloid; a dashed line, instead, 
indicates dn/vJr. „ > and therefore a local maximum of 
the critical Casimir potential with respect to D, which occurs 
within the shaded regions in Fig. [10] For a given value of A 
(with Aj < A < Ao as we shall discuss in detail further be- 
low), e.g., A = 0.60 in Fig. [TUta), the corresponding curve 
for @o shows a bifurcation at V P = X P* (A ) such that a vertical 
line drawn in Fig. [10] at a certain V P intersects this curve in 
two points @o jH and ®q s > @o,„ if *P < *F*(A), whereas it has 
no intersection for I* > *P*(A). In the former case ©o.„ and 
©0 i-s correspond to a local maximum and to a local minimum 
of the critical Casimir potential at distances Do.,, = ^ + @o.„ 
and Do s = % + ®q s , respectively, i.e., to an instable and a 
.stable levitation point for the colloid, respectively. Instead, 
for \P > *F*(A), the critical Casimir force has no zero at any 
finite value of D. We note that D = (stiction) and thus = 
always corresponds to the global minimum of the potential 
because for D — >• the critical Casimir potential is strongly 
attractive. The corresponding geometrical configuration into 
which the colloid is finally attracted by the substrate [due to 
(«i) = (— ), (b) = (— ), and X = 0, see Fig.l] is stabilized 
by the steric repulsion of the wall. We note that within the 
DA the critical Casimir potential for X — is attractive at suf- 
ficiently small distances, even if the major part of the sub- 
strate is characterized by (+) BC, i.e., even if 7^ A -c 1. 
Indeed, in this case the potential of the colloid at X = and 
close to a periodically patterned substrate can be approxi- 
mated by the one due to a single chemical lane centered at 
X = 0, which has been discussed in Sec. [V] For given col- 
loid radius R and width L\ = XP > of the attractive stripe, 
the scaling variable A = L\/(2vKD) diverges as D — > 0, so 
that the scaling function aty(A,E,®,A) which characterizes 
the potential of the lane [see Eq. ( f28b l attains the value — 1 
corresponding to the case of homogeneous, attractive (— ,— ) 
BC [see Fig. |4). Within this approximation and for D -c £± 
the critical Casimir force becomes attractive if # p ~ fig < 
which, due to Eqs. (128b . (fTTl i. and (O, yields the condition 
tBk(A,S = O,0->O,A->O)<l-2A (+ ,_ ) /(A (+ _ ) -A (+i+) ), 
i.e., cog < -0.6 in d = 4 35 and a>{ < -0.76 in d = 3 52 ; this 
occurs for A > Ao = 1.1 in d = 4, and A > Ao = 2.7 in 
d = 3, respectively [see also Fig.|4ja)]. Accordingly, at dis- 
tances D < X 2 P 2 /(ARAl) (together with D < | ± ) the criti- 
cal Casimir potential <J> p is negative and diverges to — °° for 
D — >• 0. (However, for very small values of A this would 
occur at distances of microscopic scale such that the scaling 
limit and thus the form of <J> p do no longer hold). Thus the 
bifurcation of ©0 at l P*(A) corresponds to a transition from 
(metastable) levitation at D = Dq >s for *P < ¥* (A ) to stiction 
at D = for ¥ > ¥*(A). For ¥ < ¥*(A) the metastable lev- 
itation minimum at Dq s is shielded from the global minimum 
at D = by a potential barrier the height of which vanishes 
for V F /• X P*(A) [see Fig.[TT). Experimentally, one typically 
varies the value of E,+ by changing the temperatur e 6 ' 7 ' 32 ! 53 and 
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leaves the geometry (A, P, and R) unchanged, which results 
in a change of *¥ via varying T. Thus, experimentally, the 
transition at *P* (A) corresponds to a de facto irreversible tran- 
sition from separation to stiction of the colloid as a function 
of temperature. 

Moreover, from Fig. [10] one can infer that for both d = 3 
and d = 4 there is a Ao such that, for 1 > A > Ao, K p has no 
zero for any choice of *P (i.e., there is no solution ©o) and 
the critical Casimir force is attractive at all distances. Within 
the DA, A) = A(+.-)/(A(+,-) — A(__\) [see also Eq. (l37j>l, 
which renders the values Ao = 0.80 in d = 4 s - 5 - and Ao — 0.88 
in d = 3 52 . In addition, from Fig. [10] one can infer that for 
Ao > A > Ai ~ 0.5 and *P < 1, &q s effectively does no longer 
depend on V P but solely on A. Accordingly, the distance 
Do s x <=+ at which the colloid stably levitates can be tuned 
by temperature upon approaching criticality. However, for 
A < Ai ~ 0.5, Oji diverges at ¥ = ¥ (A) < ¥*(A) such that 
for l'o(A) < ^ < l P*(A) the colloid exhibits critical Casimir 
levitation at a local minimum of the potential, whereas within 
this range of A values for ^P < ^(A) the critical Casimir 
potential has only a local (positive) maximum at Dq u ; it is 
repulsive for D > Dq m and therefore for large values of D 
(i.e., 0>1 and II ^; 1) it approaches zero from positive val- 
ues. This qualitative change in the behavior of the critical 
Casimir potential occurs at A = Ai. The value of Ai is close 
to 0.5 because the repulsive and attractive forces for (+,— ) 
and (— , — ) BC, respectively, have similar strengths but oppo- 
site signs for > 1, i.e., fc( + _)(® > 1) ~ -fc(_ _)(© > 1) 
for both d = 3 and d = 4 [see Eq. ©, where |A-'/A+| ~ 1.2 
in d = 3 7 and |A_/A + | = 1 in d = 4 35 ]. Accordingly, de- 
pending on A being larger or smaller than Ai ~ 0.5, the area 
covered by one of the two BC prevails and the resulting force 
is asymptotically (i.e., 6>1) attractive or repulsive, respec- 
tively [see the remark at the end of Sec. I VI AI and Eqs. d37] i 
and JA29H . Taking into account the slight difference in the 
strength of the asymptotic forces for (+,— ) and (— ,— ) BC 
one finds Ai = (1 — A + /A_) _1 which renders Ai = 1/2 in 
d = 4 and Aj ~ 0.545 in d = 3. The asymptotic behavior of 
the force at large distances can be inferred from the asymp- 
totic behavior of K p (X,Tl = ^©- l / 2 ,Z = 0,©> l,A->0)~ 
sf(*f?,k) ® d ~ l e~ & , which can be obtained from Eqs. ( 1351 ). 
( T33l , dA2"9T >, (|20l , and © . Accordingly, the value ^(A) at 
which ©o,.s diverges is characterized by the fact that srffi? *j 
^o (A ) , A ) ^ so that the force approaches zero from above or 
from below depending on having V P < ^(A) or V P > ^(A), 
respectively. The condition £/^¥q(X),X) = yields the fol- 
lowing implicit equation for l'o(A): 

2A^J^erf{M (n + | ) }_erf{M (n _| ) }. (4 4) 

For A <C 1 the sum on the rhs of Eq. d44l can be approxi- 
mated by the term n = alone and one finds ^(A -C 1) — 
2 3 / 2 A -1 erf (Ai), where erf - is the inverse error function, 
which yields the relations ^(A < 1) ~ 1.49/ A for d = 3 
and ^(A < 1) ~ 1.35/A for d = 4. On the other hand, in 
the marginal case one expects ^(A = Ai) = 0. However, as 
argued above, at the critical point (0 = 0) the colloid does 
not exhibit stable levitation for any geometrical configura- 
tion; this is in accordance with Fig. QJJ because for T — > T c , 
the levitation minimum of the potential moves to large D 
(Do, s = ©o,4-4+ —* °°) an d disappears at T = T c . 



In summary, as function of A there are three distinct levita- 
tion regimes: 

(i) A > A) with Ao(d = 3) ~ 0.88 and Ao(d = 4) = 4/5: 
There is no levitation and the critical Casimir force is 
attractive at all distances for any temperature. 

(ii) Ao > A > A! with A 3 (d = 3) ~ 0.545 and A 3 (d = 4) = 
1/2: Sufficiently close to T c , i.e., for »P = P/i/H^. < 
V P*(A) there is a local critical Casimir levitation mini- 
mum. Upon approaching T c its position Dq s = &q s £, + , 
with ©o,. v (4+ — > °°) finite, moves to macroscopic values 
proportional to the bulk correlation length. 

(iii) Ai > A: As in (ii) there is a local critical Casimir 
levitation minimum sufficiently close to T c , i.e., for 
*P < l F*(A). In general the onset of its appearance 
occurs further away from T c upon lowering A. Upon 
approaching T c the position Dq, s of this minimum di- 
verges at a distinct nonzero reduced temperature given 
by >Fo(A), i.e., at 4+ = P 2 /[R%(X)]: D , s - ©o,^+ 
wifh0o,*(*\YoW)-K». 

We note that, according to Figs. [6l [7] [8] and [9] we expect 
that for II < 2 and < 4 and for II < 0.5 and > 4, the 
DA does not provide a quantitatively reliable description of 
the actual behavior of K p and therefore of F p ; thus, for values 
of ¥ < 2V©0 for ©o < 4, and ^ < 0.5 % /©^ for O > 4, we 
expect quantitative discrepancies between the actual behavior 
and the one predicted by the DA shown in Fig.[l0] Nonethe- 
less our results demonstrate that the geometric arrangement of 
the chemical patterns allows one to design the normal critical 
Casimir force over a wide range. 

Figures [TTT a) and (b) show the critical Casimir potential 
Op as a function of D in d = 3 within the DA based on Monte 
Carlo simulation data for the film geometry 52 for a variety of 
specifically chosen values of the parameters P, L\, R, and |. 
The choice of these values is motivated by the typical exper- 
imental parameters which characterize recent investigations 
of the critical Casimir force acting on colloids immersed in 
binary liquid mixture a 6 i 7 i 32 ' 53 . In particular, concerning the 
colloid radius we focus on the data of Ref. [531 correspond- 
ing to R — 1.35/xm, while for the pattern we have chosen 
a periodicity P = 1/zm with A = 0.4 (i.e., L\ = 400nm and 
L 2 = 600nm) [Fig.QjIa)], or P = OAjxm with A = 0.65 (i.e., 
L\ = 260nm and L 2 = 140nm) rFig.fTTTb")!. A chemically pat- 
terned substrate with these characteristics appears to be real- 
izable with presently available preparation techniques 32,54,55 . 
[We note that 4> p as shown in Fig. [TTTa) and (b) is expected to 
describe the actual interaction potential in the scaling regime 
characterized by values of D and 4+ much larger than mi- 
croscopic length scales (such as t^ ~ 0.2nm 6,7 ) so that this 
prediction for <J> p is valid only for D,^ + > 5nm.] With this 
choice of parameters we have calculated <I> p for various val- 
ues of 4+ within an experimentally accessible range 6 ' 7 ' 32 ' 53 . 
From Figs. [TTT a) andfTTTb) one can infer that for small values 
of 4+ (corresponding to large values of VP > 1 , *(A)) the crit- 
ical Casimir potential is always attractive with a monotonic 
dependence on D [see also Fig.QJj). Upon approaching criti- 
cality, i.e., for increasing values of £, + and decreasing values 
of *P < V P* (A), a local maximum and a local minimum of the 
potential develop, so that for very small as well as for large 
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FIG. 11. Critical Casimir potential <J> p [Eq. | |32H in d = 3 of a col- 
loid of radius R = 1 .35/xm close to a periodically patterned substrate 
as a function of D and for various values of £, + for P = lfim with 
X = 0.4 in (a) and P = 0.4/im with X = 0.65 in (b) and (c). The 
values of P, X, and £, + are chosen as to be experimentally accessible 
in a colloidal suspension exhibiting critical Casimir forces ,7 ii . 
The critical Casimir potential for the colloid close to a patterned sub- 
strate may exhibit - depending on the value of §+, and, thus, on the 
temperature - a local minimum corresponding to stable levitation. 
In (c) an electrostatic potential <I> e i [Eq. {45}] is added to <l>p, which 
refers to actual experimental data . The shaded area indicates the 
ranges of the positions and the depths of the local minima of the to- 
tal potential occurring if the substrate is laterally homogeneous and 
purely attractive, i.e., for X = 1 ((—,—) BC) for the range 14nm 
< §_|_ < 75nm leading to potential depths between 0.5/^7 and 70kgT 
(indicated by the shaded arrow); for X = 1 the preferred colloid po- 
sition is dictated by the electrostatic repulsion and restricted to the 
range of 50nm to 75nm, whereas the colloid position D[). s = ©0.^+ 
due to critical Casimir levitation can be much larger and tuned by 
temperature. Moreover, whereas for X = 1 and upon approaching T c 
the minima monotonically become deeper, the levitation minima first 
deepen and move to smaller values of D followed by a decrease of 
the depth, by becoming more shallow, and moving to larger values of 
D. Reducing the range and strength of the electrostatic repulsion by 
adding salt to the solvent is expected to provide access to even more 
details of the critical Casimir levitation potential <5 p shown in (b). 



D the colloid is attracted to the patterned substrate, whereas 
within an intermediate range of values for D it is repelled 
from it [see also Fig. [TOl . Thus, the colloid stably levitates 
at a distance Dq s corresponding to a local minimum of the 
potential. The depth of this minimum ranges between a few 
kgT [Fig. [TlTa)l up to several kgT [Fig. fTTTb)!. Upon in- 
creasing §.)-, Z)o,i increases as well, i.e., the colloid position 
is shifted away from the patterned substrate with the poten- 
tial minimum becoming more shallow. In Fig. [TIT a) A = 0.4 
and we find »F(A = 0.4) ~ 4.65 and ¥ (A = 0.4) ~ 3.71 [see 
Fig.EJb)] so that for ^ < ¥ (A = 0.4), i.e., for |+ > 53.5nm 
[Fig. [TlT a)l the colloid does not exhibit stable levitation and 
the critical Casimir potential has a local maximum only. The 
levitation minimum moves to macroscopic values of D upon 
approaching the temperature corresponding to | + ~ 53.5nm. 
In Fig.QjIb) A = 0.65 and one has ¥* (A = 0.65) ~ 2.63; here 
©o,s remains finite for 'P — > in contrast to the case A < 0.545 
[Fig. [T0Tb)l. Thus, within the DA, for the case shown in 
Fig. [TTT b) stable levitation of the colloid is preserved for all 
finite values of B,+ >P 2 /[R (¥*(A =0.65)) 2 ] ~ 17nm. In this 
case upon approaching T c the levitation minimum moves to 
macroscopic values of D proportional to the bulk correlation 
length £+. 

The discussion above focuses on the position of mechani- 
cal equilibrium of the colloid, corresponding to the point at 
which the forces acting on the particle vanish and the as- 
sociated potential <I> has a local minimum 4> mm . However, 
due to the thermal fluctuations of the surrounding near-critical 
fluid at temperature T, the colloid undergoes a Brownian dif- 
fusion which allows it to explore randomly such regions in 
space where the potential <I> is typically larger than <I> m in 
for at most few IcbT. As a result, a position of mechani- 
cal equilibrium is stable against the effect of thermal fluctua- 
tions only if the potential depth of the minimum is larger than 
few IcbT. In particular, if the potential barrier <P(L\,P,0,D = 
D 0>U ,R,T) -®(Li,P,0,D = D j,R,T), which separates the 
position of the local minimum at distance D = Dq s (levita- 
tion) from the global one at D = (stiction), is not sufficiently 
large [see, e.g., the curves corresponding to £, + = 36nm in 
Fig. [TTTa) or corresponding to £, + < 18nm in Fig. [TTTb)l. a 
de facto irreversible transition from levitation to stiction may 
occur as a consequence of thermal fluctuations. 

In Fig. fTTT c) we show the resulting total potential of the 
forces acting on the colloid in the presence of an additional 
electrostatic repulsion which is experimentally practically un- 
avoidable, in order to study its effect on critical Casimir lev- 
itation. We assume that the electrostatic repulsion is laterally 
homogeneous and that it can be simply added to the critical 
Casimir potentia l 7 ' 31 ' 53 [see also Sec. lVIHl belowl. Concern- 
ing the spatial dependence of the electrostatic repulsion we 
consider the one of Ref . p3T which corresponds to a colloid of 
radius R = 1.35/xm immersed in a near-critical water-lutidine 
mixture and close to a substrate exhibiting critical adsorption 
of water or lutidine 53 : 



& d (D)/k B T = exp{-K(D-D )}, 



(45) 



where Dq = 88nm and K" -1 = llnm 53 . (Formally, 4> e i in 
Eq. J45] l is finite for D — >• 0, and thus <J> p + 4> e i is negative 
for D < 2nm and has a global minimum at D = because 
4> p -> —oo for D -» 0, However, Eq. (05]) is actually the 
asymptotic form of the electrostatic interaction which is valid 
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for distances larger than the electrostatic screening length, i.e., 
D^$> KT x . The corresponding total potential <J> p + <J> e i is there- 
fore not accurate for small values of D and is reported in 
Fig. [TTT c) for D > 50nm only.) As in Fig. [TTT b) we choose 
P = 0.4/im, A = 0.65, and experimentally accessible values 
of 4 + . Figure [TTT c) provides a realistic comparison of the 
critical Casimir potential with other forces as they typically 
occur in actual experimental systems. One can infer from 
the graph reported in Fig. [TTT c) that for this choice of pa- 
rameters the critical Casimir levitation exhibited by the col- 
loid is rather pronounced even in the presence of electrostatic 
interaction. Far from the critical point (^ + = lOnm) the in- 
teraction of the colloid with the substrate is completely domi- 
nated by electrostatic repulsion. Upon approaching criticality 
(lOnm < £, + < 35nm) a minimum in the total potential devel- 
ops and becomes deeper due to the increasing critical Casimir 
attraction working against the electrostatic repulsion. For this 
latter range of values of £,+ the local minimum of the crit- 
ical Casimir potential corresponding to levitation is located 
at distances Dq s < 60nm at which the electrostatic repulsion 
still strongly contributes to the resulting total potential [see 
Fig. [TTT c)l. Closer to the critical point (£, + > 45nm) the lev- 
itation minimum of the critical Casimir potential occurs at 
distances Dq s > lOOnm [see Fig. [TTT b)l at which the elec- 
trostatic force acting on the colloid is weak. Thus, here the 
critical Casimir effect dominates and the position of the min- 
imum of the total potential increases with increasing values 
of <^ + , which allows for measurements of the critical Casimir 
potential for distances at which the precise form of 4> e i is not 
important. Moreover, the depth of the minimum decreases 
upon approaching criticality and the minimum becomes more 
shallow. This behavior of the levitation minimum is distinct 
from the critical Casimir effect acting on a colloid close to 
a homogeneous substrate: a local minimum also occurs in 
the latter case if the critical Casimir force is purely attrac- 
tive (A = 1, (— ,— ) BC) and works against the electrostatic 
repulsion 6 - 7 , due to the competition of different forces with 
opposite sign. (We note that the critical Casimir levitation de- 
scribed above emerges from the critical Casimir force alone, 
i.e., it is a feature of a single force contribution.) However, in 
this homogeneous case the preferred colloid position D .(_,_) 
depends crucially on the form of the electrostatic interaction 
and is almost constant (50nm < Dq ,(_._) < 75nm). Moreover, 
the depths of these latter minima monotonically increase as a 
function of of E, + and become much larger than those shown 
in Fig. [TTTc) (see, e.g., Fig. 2(a) and Fig. 2(c) in Ref. |g and 
Fig. 3 in Ref.l53l). In Fig. [TTT c) this is indicated by the shaded 
area and the shaded arrow, which corresponds to the area of 
the graph within which minima of the total potential in the 
homogeneous case A = 1 occur for 14nm <<!;+< 75nm cor- 
responding to potential depths of O.SksT up to IOIcbT. On the 
other hand, the colloid position Do_ s due to critical Casimir 
levitation can be much larger, can reach values of several E, + , 
and can be tuned by temperature according to Dq, s = ®o,j<i;+. 
In conclusion, the examples presented in Fig. [TT] strongly sug- 
gest that the critical Casimir levitation of a colloid close to a 
patterned substrate is experimentally accessible. 

By patterning the substrate, one introduces an additional 
(lateral) length scale into the system, which, according to 
our results presented above, can finally lead to stable levita- 
tion. Introducing an additional length scale along the nor- 



mal direction by stacking different materials on top of each 
other may lead to levitation due to quantum-electrodynamic 
Casimir forces 56 . The behavior of the stable levitation dis- 
tance shows a bifurcation and irreversible transitions from 
separation to stiction 56 similarly to the ones described above 
[see Fig. [TO]. In that context great importance has been given 
to the temperature dependence of the position Dq, s of stable 
quantum Casimir levitation 56 , which is quantified by the value 
of -jjDq s . In the critical Casimir case presented here, for an 
estimate of -jyDq s we pick as an example the stable levitation 
positions for | + = 18nm and <i; + = 60nm as reported in Fig.Qj] 
(a different choice would lead to similar results). The results 
reported in Fig. QT| correspond to the experimentally relevant 
water-lutidine mixture with £j" = 0.2nm and T c ~ 307K 6 ' 7,53 . 



Therefore, according to 4+/<jj 



Tc)/T c 



" v , the dif- 

= 18nm to 
-l 



ference in temperature required to move from ^ 
£,+ = 60nm is AT ~ 0.2K. Thus we find ^D , s ~ 560nm K 
for the average temperature dependence of critical Casimir 
levitation [Fig.[TTTb)l, and ^D 0tS ~ 230nm KT 1 by addition- 
ally taking electrostatics into account [Fig. [TTT c)l. We note 
that in the present critical case -jL Dq s can become arbitrar- 
ily large at temperatures corresponding to the transition from 
separation to stiction and the emergence of the local mini- 
mum and the local maximum of the critical Casimir potential 
[see Fig.fTOland the curves for <i; + = 34nm and ^ + = 36nm in 
Fig. [TTT a)l. This shows that the critical Casimir levitation is 
strongly temperature dependent, even near room temperature, 
with the variation of stable separation ;jf A),.« being two orders 
of magnitude larger than the one predicted for the quantum- 
electrodynamic Casimir effect in Ref. \5($. In general the col- 
loid will not only be exposed to the critical Casimir force and 
to an electrostatic force but also to gravity and to laser tweez- 
ers, which generate a linearly increasing potential contribu- 
tion. This attractive contribution tends to reduce the potential 
barriers shown in Fig. [TT] and can eliminate small barriers al- 
together. Thus these external forces can be used to switch 
levitation on and off (compare a similar discussion related to 
the quantum-electrodynamic Casimir levitation in Ref. do). 



VII. CYLINDER 

Currently, there is an increasing experimental interest in 
elongated colloidal particles which have a typical diameter 
of up to several 100 nm and a much larger length (see, e.g., 
Refs.l30landl57l and references therein). These types of col- 
loids resemble cylinders rather than spheres. The description 
of their behavior in confined critical solvents calls for a natu- 
ral extension of the studies presented in Secs. HIH4Vll Hence, 
in the present section we consider the case of a 3d cylinder 
with (— ) BC which is adjacent and parallel aligned to a pe- 
riodically chemically patterned substrate consisting of alter- 
nating (— ) and (+) stripes as the ones discussed in Sec. I VII 
Accordingly, the axis of rotational invariance of the cylinder 
is perpendicular to both the x direction [Fig. [JJ and the di- 
rection normal to the substrate, and it is parallel to the direc- 
tion of spatial translational invariance of the chemical stripes 
forming the pattern. As compared with the case of the sphere 
the analysis for the cylinder is technically simpler because the 
system as a whole is invariant along all directions but two, 
the lateral one, x, and the one normal to the substrate. (For 
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the sphere its finite extension in the second lateral direction, 
which is normal to the jc-axis, matters and thus leads to a ba- 
sically three-dimensional problem. Accordingly, here we do 
not consider short cylinders, for which this finite length mat- 
ters, too.) This reduction of the number of relevant dimen- 
sions allows us to perform numerical calculations of adequate 
precision for a range of various pattern geometries which is 
wider than in the case of the sphere. (Here, we do not con- 
sider a cylindrical colloid which is not perfectly aligned with 
the pattern and which would, therefore, experience a critical 
Casimir torque 30 .) Even though the expressions derived in 
Appendix iDl can be used to study the case of a cylinder hav- 
ing its axis laterally displaced by an arbitrary amount X from 
the chemical step, our numerical calculations for the case of a 
chemical stripe address only the case X = 0. This corresponds 
to a lateral position of the symmetry axis of the cylinder which 
coincides with the center of an attractive (— ) stripe. 

In Appendix [D] we briefly derive the scaling behavior of 
the normal critical Casimir force acting on the cylinder and 
compare it with the case of a sphere. Then, we adapt the 
Derjaguin approximation appropriate for the geometry of the 
cylinder. On this basis, we have calculated the scaling func- 
tion of the normal critical Casimir force acting on the cylinder 
in d = 3 and d = 4 on the basis of the Monte Carlo simulation 
data for the film geometry 52 and of the analytic MFT expres- 
sion for the critical Casimir force for the film geometry^ 5 -, re- 
spectively. In addition, within the same approach as the one 
of Sec. Ill Bl we have calculated numerically the MFT scaling 
functions corresponding to A ^ 0, in order to assess the per- 
formance of the DA. 

Here we focus on the comparison between the DA appro- 
priate for the cylinder and the full numerical MFT data for the 
scaling function Kp y (A,n,E — 0,©,A) which characterizes 
the normal critical Casimir force in the presence of a peri- 
odically patterned substrate; A, II, S, 0, and A are defined 
as in the case of the sphere [see Sec. IVII and Appendix ID1 . 
Figure [12] shows the scaling function of the normal critical 
Casimir force acting on a cylinder as a function of © as ob- 
tained from the DA (A — > 0) in d = 4 and from the full nu- 
merical MFT calculations for A = 1/3. Besides the quanti- 
tative differences in the scaling function as a function of ©, 
the qualitative features of the behavior of the force acting on 
a cylinder, which is reported in Fig. Q~2] for various values of 
A, are similar to the ones for the sphere [compare Fig. [£)■ 
For n = 1.92 [Fig. []~2 a)] the DA describes the actual behav- 
ior of the critical Casimir force rather well, in particular for 
© > 2, even for most values of A. As in Fig. [9] for a certain 
range of values of A the normal critical Casimir force changes 
sign at ©q Y (IT, A, 3 = 0, A). On the other hand for small pe- 
riodicities (II = 0.29 in Fig. fl27b)) the DA in d = 4 fails to 
describe quantitatively the actual behavior of the force as ob- 
tained from the full numerical MFT calculations. These strong 
deviations from the DA rFig. [T2T b)l indicate the relevance of 
effects caused by the actual non-additivity of critical Casimir 
forces. 

For A > 0.6 the scaling function K^ y of the normal crit- 
ical Casimir force obtained numerically and represented by 
symbols in Fig. [T2T b) is very close (much closer than within 
the DA) to the one corresponding to the homogeneous case 
with ( — , — ) BC (corresponding to A = 1) and does not show 
a change of sign. This means that, even if the substrate is 
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FIG. 12. Normalized scaling function A*p Y [see Appendix IDl includ- 
ing expressions for K^ _, (0, 0)] of the normal critical Casimir force 
acting on a cylindrical colloid close to and parallel to a periodically 
patterned substrate. The cylinder axis is aligned with the striped pat- 
tern and positioned above the center of a (— ) stripe which has the 
same adsorption preference as the cylinder (analogous to Fig.[9]for a 
spherical colloid). In (a) for II = 1.92 the appropriate DA describes 
the actual MFT data rather well, and for 0.3 < X < 0.7 there is a 
change of sign of the force. In (b), instead, apart from the limiting 
homogeneous cases A = and X — 1, for n = 0.29 the DA fails to 
describe quantitatively the actual behavior [see the main text]. In (c) 
Kp y is shown for d = 3 within the DA based on the Monte Carlo 
simulation data for the film geometry 52 for the two cases IT = 0.29 
and II = 1.92. We expect that also in d = 3 the DA for n = 0.29 is 
not quantitatively reliable. 



not homogeneous but chemically patterned - but such that the 
larger part of the surface still corresponds to (— ) BC, i.e., 
A > 0.5 - the resulting critical Casimir force acting on the 
colloid with (— ) BC resembles the behavior for laterally ho- 
mogeneous (— ,— ) BC. This can be understood in terms of 
the fixed point Hamiltonian in Eq. © which penalizes spa- 
tial variations of the order parameter at short scales. Thus 
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the system tries to smooth out spatial inhomogeneities of the 
order parameter profile, biased by the preference of the col- 
loidal particle. If the pattern is very finely structured, i.e., 
II = [L\ +L2)/ \/RD -C 1, regions with a positive order pa- 
rameter close to the narrow (+) stripes (A ~ 1 _ , i.e.,L2 ^^l) 
extent only very little into the direction normal to the substrate 
and the resulting order parameter profile at a distance from 
the substrate remains negative only 58 , so that the force resem- 
bles the one corresponding to the homogeneous case. (Note 
that within the DA, the corresponding order parameter profile 
would simply consist of a patchwork of the order parameter 
profiles corresponding to the film geometry, with no smooth- 
ing taking place at the edges of the various spatial regions.) 
Similarly, but in a weaker manner due to the opposite order 
parameter preference at the colloid, the curves in Fig. [T2T b) 
for A < 0.5 approach the corresponding homogeneous one 
for the case (+,— ) (i.e., A = 0). Thus, the fact that both in 
Fig. fTZl' a) and Fig. [T2"l" b) the curves for X = 1/5 are less close 
to their limiting ones for X = than the curves for A = 4/5 
are close to the ones for A = 1 - although the portions of the 
minority part of the surface are the same - is due to the fact 
that an order parameter profile with (+, — ) boundary condi- 
tions is energetically less preferred than the one with (— ,— ) 
boundary conditions because in the (+,— ) case an interface 
emerges between the two phases. For broad stripes, i.e., in 
contrast to the case n — > 0, the energy costs for a similar be- 
havior are seemingly larger: the full numerical MFT data for 
X = 1 /5 and A = 4/5 are less close to the corresponding lim- 
iting homogeneous cases A = and A = 1, respectively, for 
n = 1.92 than for n = 0.29. 
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Figure fT2T c) shows the scaling function K^ y of the normal 
critical Casimir force for d = 3 within the DA as obtained by 
using Monte Carlo simulation data for the film geometry 52 . 
One can infer from Fig. fT2l c) that the qualitative features 
of the MFT scaling function as described above, such as the 
change of sign, are carried over to d = 3. 

As discussed in the previous section, the vanishing of the 
normal critical Casimir force corresponds to a stable levita- 
tion of the colloid at a distance Do from the substrate only if 
doFp yi I d=d < 0. Within the DA and at the laterally stable 
position E — the sign of dj)F ? cy is given by Eq. (1431 with 



^eyl 



The behavior of © C J as a function of 



Kp replaced by K^ 
*P and the demarcation of the regions where levitation is sta- 
ble against perturbations of D is shown in Fig. [13] where the 
solid and the dashed lines correspond to stable and unstable 
levitation, respectively. The behavior for the normal critical 
Casimir force acting on the cylinder is qualitatively similar to 
the one for the sphere shown in Fig. [TOj Analogously to the 
case of a sphere discussed in Sec. I VI CI no stable levitation 
is found at T = T C or for A > Xq — Ar + \/(Ar + \ —Ar\), 
where Ao = 0.80 in d = 4 and Ao ^ 0.88 in d = 3. On the 
other hand, for > 0, and A < Ao, it is always possible to 
find values of P and R such that stable levitation of the cylin- 
der occurs at a certain distance from the substrate. The val- 
ues of Ai below which one has a finite value x ¥q(X) at which 
©0 diverges remain the same as for the case of a sphere, i.e., 
X\ (d = 4) = 1/2 and X\ (d = 3) ~ 0.545; also the correspond- 
ing values of x ¥q(X) remain the same [see Eq. (l44l l. 
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FIG. 13. Values of the scaling variable Off at which the normal crit- 
ical Casimir force A^ y acting on a cylinder close to a periodically 
patterned substrate vanishes as a function of *P = P/ \/R%+ [com- 
pare Fig. [TO] for the case of a sphere] within the DA. The region 
indicated by solid lines corresponds to the one in which the levi- 
tation of the cylinder at a height D = Do = &q£,+ is stable against 
small perturbations of D, wheres in the shaded region indicated by 
dashed lines there is no such stable levitation although the normal 
critical Casimir force acting on the colloid vanishes. For X > Ao 
with l (d = 4)= 4/5 and X Q (d = 3) 



0.88, 0Q y ceases to exist, 



i.e., A" p does not exhibit a zero. For A < X\ with X\ (d = 4) = 1/2 
and X x (d = 3) ~ 0.545, ©jf^ \ *o W) diverges. (The values for 
x Vq(X) are indicated by upward arrows.) For any X < Aq, 0q Y exists 
for *P < I** (A). We expect the DA to be quantitatively reliable only 
for 'P/v 7 ©!) > 2 for O < 4 and for ^/V®o £ 0.5 for G > 4. 



VIM. SUMMARY AND CONCLUSIONS 

We have investigated the universal properties of the nor- 
mal and lateral critical Casimir forces acting on a spherical or 
cylindrical colloidal particle close to a chemically structured 
substrate with laterally varying adsorption preferences for the 
species of a (near) critical classical binary liquid mixture (at 
its critical composition) in which the colloid is immersed. 
Within the Derjaguin approximation (DA) [see Fig. [14) in spa- 
tial dimensions d = 3 and d = 4 we have derived analytic ex- 
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pressions for the corresponding universal scaling functions of 
the forces and the potentials for general fixed-point boundary 
conditions (BC) in terms of the scaling function of the criti- 
cal Casimir force acting on two parallel, homogeneous plates. 
These expressions are given explicitly analytically at the bulk 
critical point T = T c and - for symmetry breaking boundary 
conditions - far away from the critical point. These relations 
enable one to obtain predictions for actual three-dimensional 
systems with a sphere-inhomogeneous plate geometry (for 
which currently computations are not possible) based on the 
scaling function for the parallel homogeneous plate geometry, 
for which, e.g., Monte Carlo simulation data in d = 3 are avail- 
able. Moreover, results within mean-field theory (MFT, cor- 
responding to d = 4) and symmetry-breaking boundary con- 
ditions rSec. Ill Bl have been obtained fully numerically and 
have been compared with the approximate results of the DA, 
which allows us to explore the limits of validity of the latter. 
We have studied several relevant situations [see Fig. Q~) and 
our main findings are the following: 

1. First, we have studied a spherical colloid immersed in 
a binary liquid mixture close to a chemically homoge- 
neous substrate which has, compared to the colloid, the 
same (— ) or a different (+) adsorption preference for 
one of the species of the mixture [Sec. [III). Close to 
the bulk critical point at T = T c the critical Casimir 
force induced by the confinement of the order param- 
eter (e.g., the concentration difference in a binary liquid 
mixture) can be described in terms of universal scaling 
functions depending on the surface-to-surface distance 
D of the colloid from the substrate scaled by the bulk 
correlation length, © = sign((T — T c )/T c )D/£,±, and its 
ratio with the radius of the colloid, A = D/R [Eqs. (0 
and ©]. The scaling functions obtained within the DA 
[Eqs. ( [Tol l and ( fTTT ll are valid for A — > 0. From the com- 
parison with the full numerical MFT results [Fig. [2] we 
find that in d = 4 the DA describes the actual behavior 
quite well for A < 0.4. Based on Monte Carlo simula- 
tion data for the scaling function of the critical Casimir 
force between parallel, homogeneous plates and within 
the DA we have obtained also the scaling function for 
the critical Casimir force on a spherical colloid close to 
a homogeneous substrate in d = 3 [Fig. [2]. 

2. The basic building block of a chemically patterned sub- 
strate is a chemical step, which we have studied in 
Sec. |IV] Due to the broken translational invariance in 
one lateral direction (x) the critical Casimir forces and 
potentials acquire a dependence on the additional scal- 
ing variable S = X/^/RD, which corresponds to the lat- 
eral distance X of the center of the spherical colloid 
from the position of the chemical step along the plane 
[Eqs. (O, (TBI , and (l2Tb l. Due to the different bound- 
ary conditions on both sides of the chemical step a lat- 
eral critical Casimir force emerges, which leads to a 
laterally varying potential for the colloid. In the limit 
A — » both the scaling function for the potential and 
for the lateral critical Casimir force as obtained within 
the DA are in agreement with the full numerical data 
[Fig. [3). We have derived the corresponding scaling 
functions within the DA also in d = 3 by using Monte 
Carlo data for the parallel plate geometry [Fig. [3)- The 



preceding results have been partly presented in Ref . |3l| 
as well as their suitable comparison with corresponding 
experimental results 32 , which revealed that the critical 
Casimir effect is rather sensitive to the geometrical de- 
tails of the substrate patterns. 

3. Section [V] deals with the critical Casimir forces and 
the corresponding potential acting on a spherical col- 
loid in front of a single chemical lane of width 2L, 
which additionally depends on a fourth scaling vari- 
able A = L/\/RD [Eqs. d25j and ([26}]. It turns out 
that within the DA the scaling functions for the critical 
Casimir force and the critical Casimir potential across a 
chemical lane can be expressed in terms of the ones for 
the chemical step [Eqs. ( |29l and (f30T>1 - For large values 
of A the resulting potential can be described as a suit- 
able superposition of chemical steps, whereas for A < 3 
one has explicitly to account for the finite width of the 
chemical stripe [Fig. 0J. Comparing the results of the 
DA with the ones obtained by a full numerical analy- 
sis, one finds that the DA describes the actual behavior 
quite well for A < 0.4, even for small A. Seemingly, 
in this respect, the nonlinearities inherent in the criti- 
cal Casimir effect and edge effects do not considerably 
affect the resulting scaling functions [Fig. [5). 

4. On the basis of the results of Sec. [V] in Sec. IVII we 
have studied the universal scaling functions of the crit- 
ical Casimir force and the corresponding potential for 
a sphere opposite to a periodically patterned substrate 
with laterally alternating chemical stripes of different 
adsorption preferences [Sec. lVII . These scaling func- 
tions [Eqs. ( |3T| i and ( |32l l depend, besides the scaling 
variables 0, A, and S, on two additional scaling vari- 
ables n = P/^/RD and A = L\/P, which correspond to 
the period P = L\ +L2 of the pattern and to the width 
L\ < P of the stripes with the same adsorption prefer- 
ence as the colloid. The scaling function for the normal 
critical Casimir force obtained within the DA can be ex- 
pressed in terms of the one for the chemical step and de- 
scribes the actual behavior well for II > 2 [Eq. d35l l and 
Figs. [3a) andUta)]. However, for II ->• [Eq. 03] 
the DA fails to capture quantitatively the numerically 
obtained behavior within MFT, reflecting the impor- 
tance of nonlinearities and edge effects in this context, 
which are not accounted for by the DA [Figs. [6] Ul 3 ) 
and gib)]. The failure of the DA in the limit II -> 
can be traced back to the fact that for the film geometry 
of a patterned wall next to a laterally homogeneous fiat 
wall, additivity of the critical Casimir forces does not 
hold [Fig. [8). 

5. The MFT scaling function of the normal critical 
Casimir force acting on a colloid close to a periodi- 
cally patterned substrate shows a remarkable behavior 
as a function of = D/£, + . Within a certain range of 
values of II and X the critical Casimir force vanishes 
at ©o corresponding to a distance D — Dq between the 
colloid and the substrate. We have analyzed the sign 
of the derivative of the critical Casimir force with re- 
spect to D at Do, which is negative if for D < Dq = Dq s 
the colloid is repelled from the substrate whereas for 
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D> D = D 0s it is attracted to the substrate [Fig.fTUl. 
This means that in the absence of other forces the col- 
loid can levitate above the substrate at a stable distance 
which can be tuned by temperature. Stable levitation 
points are found also in d = 3, within the DA and on the 
basis of the Monte Carlo data for the parallel plate ge- 
ometry [Figs.[7Jb),[9jc), and [T0f b)l. Our analysis shows 
that at the critical point T = T c levitation is not possi- 
ble, whereas off criticality a geometrical configuration 
leading to stable levitation can always be found. For 
fixed geometrical parameters, the critical Casimir po- 
tential as a function of D changes from a monotonic be- 
havior to a non-monotonic one upon approaching criti- 
cality; a local maximum and a local minimum, the latter 
corresponding to stable levitation, occur [Fig.[TTJa) and 
(b)]. Experimentally, this corresponds to a de facto ir- 
reversible transition from separation to stiction of a col- 
loid and a patterned substrate. The depths of these po- 
tential minima can be up to several IcbT so that the levi- 
tation is stable against Brownian motion of the colloid. 
The critical Casimir levitation can be rather pronounced 
and robust even in the presence of electrostatic interac- 
tions [Fig. [TTT c)l. The levitation height is proportional 
to the bulk correlation length and thus can be tuned 
by varying temperature. Depending on the geometric 
parameter A we have identified two distinct types of 
temperature dependences of the levitation height Do.. s . 
In both cases it exhibits a high temperature sensitivity 
-jjDq s which, for realistic examples at room tempera- 
ture, is of the order of several lOOnm K _1 . These re- 
sults show that the periodic patterning of the substrate 
enables one to design critical Casimir forces over a wide 
range of properties. 

6. This behavior is also observed for a cylindrical colloid 
which lies parallel to the substrate such that its axis is 
aligned with the translationally invariant direction of the 
stripes [Sec. lVlll and AppendixlB"! . The main features of 
the scaling function for the corresponding normal criti- 
cal Casimir force are similar to the ones for the spheri- 
cal colloid: the DA describes well the actual behavior as 
obtained from full numerical MFT calculations for large 
values of n, but fails quantitatively for II < 2 [Fig. [T2l . 
The numerical studies for II — > indicate that a sub- 
strate with a very fine pattern, dominated by one of the 
two BC as far as the corresponding covered area is con- 
cerned, leads to a normal critical Casimir force which 
resembles the one for a homogeneous substrate char- 
acterized by the dominating BC [Fig.[T2Tb)l. Based on 
Monte Carlo data for the parallel plate geometry we cal- 
culated within the DA the critical Casimir force acting 
on a cylinder in d — 3 [Fig. [T2l' c)l. Above a chemically 
patterned substrate, also for a cylinder stable levitation 
is possible for a wide range of parameters [Fig. [13). 

Typically, in experiments with a colloidal suspension one 
has to consider also other forces, such as electrostatics, grav- 
itation, and van der Waals forces which act on the colloidal 
particles in addition to the critical Casimir forces. The to- 
tal force is approximately the sum of these contributions^! , 
[see Fig. [TTTc)l. Upon approaching the critical point in 
the phase diagram, experiment s 6 ' 7 ' 32 and theory (see, e.g., 



Refs. [J [3j], and 1331) highlight the importance and the rele- 
vance of the critical Casimir effect in comparison with these 
other forces. 

The lateral critical Casimir forces occurring for patterned 
substrates as discussed here are highly sensitive to the details 
of the geometry of the pattern. A detailed comparison with 
available experimental data 32 has to take this into account 31 . 
This sensitivity even allows for an independent determina- 
tion of the geometry of a chemically structured substrate by 
means of the critical Casimir effect. This is useful in cases 
in which it is difficult to infer the geometry of the chemical 
pattern directly 31,32 . Concerning the comparison with experi- 
ments for chemically structured substrates, the theoretical pre- 
dictions for the critical Casimir force are in agreement with 
the presently available dat a 31 ' 32 , for which the description in 
terms of independent chemical steps [Sec. llVII turns out to be 
sufficient 3 1 . In order to test our specific predictions obtained 
for narrow single chemical lanes and for periodic chemical 
stripes, structures on the nanometer scale are needed. Prelim- 
inary experimental data in this direction are encouraging 54,55 . 

In view of present basic research efforts and potential ap- 
plications, it is important to study the effect of weak critical 
adsorption of the fluid at the confining surfaces, correspond- 
ing to finite surface fields. Such weak surface fields can be 
realized by applying suitable surface chemistry and they in- 
fluence the resulting behavior of the critical Casimir effect 
strongl y 15 ' 53 . Another approach to create an effective reduc- 
tion of the surface adsorption is to create fine periodic chem- 
ical patterns with different (strong) adsorption preferences as 
discussed here. However, our results for II — >• [Figs. |9jb) 
and [T2T b)l show that a fine patterning of the substrate with 
alternating boundary conditions does not necessarily lead to 
an effective reduction of the surface adsorption at short dis- 
tances because in this range the critical Casimir force for a 
inhomogeneous adsorption preference resembles the one for 
a homogeneous substrate corresponding to strong adsorption. 
On the other hand, at large distances a periodically patterned 
substrate does lead to an effective BC corresponding to a weak 
adsorption preference, and for X = 1/2 the surface fields even 
cancel out, leading to an effective BC resembling the so-called 
ordinary BC 33 . This offers the interesting perspective to study, 
at least asymptotically, critical Casimir forces with Dirich- 
let BC by using classical fluids instead of superfluid quantum 

fluids 10,17-18-M . 

A patterning on the molecular scale is not captured by the 
continuous approach pursued here, which gives the universal 
features of the critical Casimir effect. Nonetheless, a molecu- 
lar patterning of the substrates may provide another means for 
an effective reduction of the adsorption of the corresponding 
fluid at the surface. However, on a molecular scale the pat- 
terning is more likely to lead to randomly distributed surface 
fields which opens a new challenge in the context of critical 
Casimir forces. 
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Appendix A: Derjaguin approximation for a chemical step 

In this appendix we first calculate within the DA the nor- 
mal critical Casimir force F S (X,D,R,T) [Eq. (fT2l l acting 
on a spherical colloid of radius R facing a chemical step 
by using the DA. (We cannot directly calculate the lateral 
critical Casimir force F s "(X,D,R,T) within the DA because 
for two parallel homogeneous plates such a force vanishes.) 
In a second step we derive the critical Casimir potential 
<P S (X,D,RJ) = f£dzF s (X,z,R,T) by integrating this re- 
sult for the normal critical Casimir force. In a third step the 
lateral critical Casimir force is obtained as Ff(X,D,R,T) — 

-d x ^(X,D,R,T) = -j;dzd x F s (X, Z: R,T)[s^Sec.MM- 
In the spirit of the DA, the surface of the spherical col- 
loid with (b) BC is thought of as being made of a pile of 
(infinitely thin) rings parallel to the opposing substrate and 
with an area dS(p) = 2npdp, where p is the radius of the 
ring. Each of these rings is partly facing (in normal direc- 
tion) the surface with (a<) BC, with an extension dS < (p), and 
partly facing the surface with (a>) BC on the other side of the 
chemical step [Fig.[l4), with an extension dS > (p), such that 
dS(p) — dS<(p) +dS > (p). For an assigned p, dS^(p) de- 
pend, inter alia, on the lateral position X of the colloid. Using 
the assumption of additivity of the forces underlying the DA 
we suppose that the contribution d,F s (p) of the ring to the to- 
tal critical Casimir force F s is given by the sum of the critical 
Casimir forces which would act, in a film, on portions of ar- 
eas dS< and d5> in the presence of (a < ,b) and (a > ,b) BC, 
respectively. According to Eq. (HJ this leads to the following 
expression for the force acting on a single ring: 



dS>(fi) 



dF s (p) _ dS < (p) 

L d {p) {a « b) 

, dS>(p) 



k B T 



U>(p) 



(sign(f)L(p)/^ ± ) 

, &) (sign(0L(p)/<^ ± ), (Al) 



where L(p) is the substrate-ring distance [Fig. fT4ll as given 
in Eq. (J7J, and kr a> y\ are the scaling functions of the critical 
Casimir force in the film geometry with («>,/?) and (a<,£>) 
BC, respectively [see Eq. (Q}]. This assumption neglects all 
edge effects along the boundary between the areas d5> (p ) and 
dS < (p), which might actually be relevant in view of the spa- 
tial variation of the order parameter profile. It is therefore im- 
portant to test the validity of this assumption at least in some 
relevant cases. This is carried out in Sec. [IV] for d = 4, i.e., 
within MFT 

Without loss of generality in the following we assume X > 
0, i.e., that the normal projection of the center of the sphere 
falls on the part of the substrate with (o>) BC [Figs. Q]and 
[Pfl . The results for X < are obtained by exchanging in the 
formulas below a< ^- a > and X <H- —X. Taking into account 
that dS(p) = dS < (p) +d5 > (p) one can rewrite Eq. ( IAU as 



dF s (p) _dS(p). 
k B T L d {p) 



jfc) (sign(f)L(p)/£±) 



dS<(p) 
U'{p) 



M(sign(r)L(p)/£t), (A2) 



where Ak(&) — kt a<l b) (®) ~ ^(«>,fe) (©)■ Summing up all force 
contributions from the rings of different radii p, one finds for 




FIG. 14. Sketch concerning the Derjaguin approximation for the crit- 
ical Casimir force acting on a colloid opposite to a chemical step. 
The critical Casimir force is subdivided into contributions from rings 
parallel to the substrate. The projection of the area dS(p) of a ring 
onto the substrate is separated into the areal contributions d5< and 
dS> which emerge as the intersection of the projected ring with the 
half-planes carrying (a<) and (a>) BC, respectively [see the main 
text]. The sphere has a surface-to-surface distance D from the sub- 
strate and its center has a lateral distance X from the chemical step. 



the total normal force acting on the sphere 

F s {X,D,R,T)=F {a>ib) (D,R,T)+AF(X,D,R,T), (A3) 

where Fi a>h \ is the force acting on a sphere close to a homo- 
geneous substrate with (a>) BC and is given by Eq. (O or by 
Eqs. <(5j and (HOI , This term does not contribute to the lat- 
eral critical Casimir force experienced by the colloid near the 
chemical step, because it does not depend on the lateral co- 
ordinate of the colloid. The second term AF in Eq. ( ]A3t cor- 
responds to the integration of the force differences Ak in the 
region of overlap between the projection of the sphere onto 
the substrate plane and that part of the substrate with (a<) 
BC. For each ring this area is given by [see Fig. [14) 



dS < (p) = 



This leads to 



0, 



P<X, 



2arccos(X/p)pdp, X <p<R. 



(A4) 



AF(X,D,R,T) 
kWf 



= 2 



dp p arccos — 
P 



X\ M(sign(r)L(p)/k) 



Li{p) 



(A5) 
In the spirit of the DA, the radius of the sphere is taken to 
be large compared to its distance to the substrate, i.e., A = 
D/R <§; 1, and the contributions from the rings closest to the 
substrate dominate. Therefore, it is well justified and in accor- 
dance with the DA to assume X/R <C 1 because the contribu- 
tions of rings with large radii do not change the behavior of the 
force in the Derjaguin limit. Within these two limits we can 
use the parabolic approximation for the distance of the rings 
to the substrate [Eq. Ql], L(p) ~ Da, with a = 1 + p 2 /2RD. 
Changing the integration variable in Eq. dA5b we directly find 



R 



AF(X,D,R, T) = kaT—^AKiZ, 0, A) 



(A6) 
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where AK is a universal scaling function given by f£ dlF s (X,l,R,T) can be separated analogously to Eq. ( lA31 l. 

i.e., 

AK(Z,®,A^0)=2 f dacT d arccos( —^ — )AJk(a©). 

J \V^^)J 9,(X t D,R,T)=* ia>Jb) (p,R,T)+A9(X t R,D,T) (A8) 

l+S 2 /2 v ' 

(A7) 
Note that the relevant scaling variable E = X j\jRD can with 
take on arbitrary values, irrespective of the two assumptions 

D/7? < 1 and X/R < 1. From Eq. (A7]i one finds with " „ 

Eqs. dTOl and O directly the expression for the scaling func- A®(X,R,D, T)= dlAF(X,l,R, T) =: k B T —j^ Ai?(S,0, A) 
tion V(a<|a > ,6) given in Eq. (16). d 

The critical Casimir potential <P S (X,D,R,T) = (A9) 



Using Eq. ( 1A71 >, the scaling function A# is given by 



1 f . 1 



A#(S,0,A)=2ydy— y da^arccosl _ 1 Ak{ya&). (A10) 

By changing the integration variable a*-> z:= 2y(a — 1)/S 2 followed by y i->- v := y + E 2 z/2 one obtains 

A$(Z,@,A) = Z 2 [ dz [ dv4arccos(l/Vz)A/t(v>0). (All) 

VI Jl+zS 2 /2 y 

After changing the order of integration 

dz / dv = / dv / dz, (A12) 

1 il+z» 2 /2 il+S 2 /2 J\ 



and using the primitivi 
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dzarccos(l/ v / i) = zarccos(l/ v / z) — y/z—l+c, (A13) 



A#(E,0,A) = ^- / ds ^ ; „ 2 y [yarccos(^ 1/2 )-V^TJAj!:(0[l+sV2]). (A14) 



one obtains after a final change of variables v !-► w := 2(v— 1)/S 

Til " J (1+E 2 V2) 

From Eq. ( IA14t together with Eq. ( fTTT i one obtains the final expression for the scaling function of the critical Casimir potential 
as given in Eq. (TT~8T >. 

1 . Bulk critical point: = 

In order to calculate the critical Casimir force acting on the colloid at the bulk critical point one inserts Eq. (ffjl into Eq. ( IA7t 
and obtains 



AAT(E,0 = O,A) = 2(A {a<>b) -A {a>tb) )J daa- d KCCOs\-j=^=) (A15) 



=:3 2 (A (fl<>fe) -A (fl>ife) )/ d (S 2 /2), 

I 

where A( a m = hub) (0) t see Eq. ©], and with the substitution For I t i{a) the recursion relation 



a ^z = Z/y / 2(a - 1) for d > 1, 



I d+1 {a) = \a l - d ^-{cfil d {a)} {All) 

a da 



-1 T 2rf-3 



Pl -AU — 3 

d ^ ' J (z 2 +a) d holds, so that I4 and 73 can be expressed in terms of h- Per- 
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Ha) = - 



forming the integration we find 63 

1- 

and therefore with Eq. ( IA17I ) 
h{a) 



,1/2 



% 

4a 



1 



(l+a)V2 



\a l l 2 + c?l 2 



and 



7 4 (fl) 



6« 



1- 



(1+fl) 3 / 2 

ffl 1 /2 + 5 fl 3/2 + fl 5/2 
(1+0)5/2 



(A18) 



(A19) 



(A20) 



Thus, from Eqs. JA16t , ( |A19b , and dA201 > together with the 
expression for Kr a> m(0,0) = 2nA/ a>jb \/(d— 1) [Sec. IIII Al 
and Eq. (fT3l . one finds the expression for the scaling function 
W(a < \a >< b) given in Eq. ( f+Tb - The critical Casimir potential at 



© = for d = 3 and 4 can be found from Eq. ( lAlOt together 
with Eq. ( IA16b : 

A#(S,0 = 0,A) = S 2 (A (fl<>6) - A (fl>)6) ) j^dyy-% (§) , 

(A21) 
and from a change of variable y4a = E 2 /(2y) one finds 

2 <*-2 |-S 2 /2 

A^(E,0,A) = ^2 Zr5 (A ( „ < , i) -A ( „ >i)) )y daa d - 2 I d (a). 

(A22) 
Using Eq. dA171 > and the limiting behavior I c [(a — > 0) = 
7r/(2(,i-l)fl),wefind 



A0(E,O,A) 



2(J-1) ( A (»<^)~ A (»>- fc )) /t/ - 1 ( :5:2 / 2 )- 

(A23) 
From Eqs. ( TA18b, dA19l ), and ( TA231 together with i? (0,0) as 
given in Sec. IIII Al one obtains Eq. dT9b for the scaling function 
of the critical Casimir potential at T c . 



2. Far from criticality: © > 1 

Far from the critical point, i.e., for © 3> 1, and for symmetry breaking boundary conditions (o<) = (+), (a>) = (— ), and 
(b) = (—) Eq. (f3]l holds and the integrals in Eqs. (IA71 > and (!A14b can be calculated analytically. For 3> 1 Eq. dA7b turns into 



A/r(S,0> l,A)=2(A_-A + )0 c/ / daarccos. 



-a& 



Substituting a i-> j3 = 2(a — 2)/E 2 one has 



A^(S,0>l,A)=S 2 (A_-A + )0"e- w / dj3arccos(j3"> 



U„-E 2 00/2 



Integrating by parts leads to 



By using the relation 64 



AK(Z 1 ®^>l,A) = {A--A+)® d - 1 e- & dj3 



-Z 2 0j3/2 



i PVP^ 



dj3- 



= e - fl2 ^ = Trerfc(a), 



(A24) 



(A25) 



(A26) 



(A27) 



where a > and erfc(a) = 1 — erf (a) = 2n: -1 / 2 f™dt exp(— t 2 ) is the complementary error function, we finally arrive at 

A^(E,0>l,A) = ^(A_-A + )0 rf - 1 e- erfc(E\/072). (A28) 

The scaling function K^—\ for © 3> 1 in the homogeneous case [Sec.|III) is given by 7 

%_)(0» l,A^O)=2^A ± J - 1 e-® (A29) 



and from Eqs. (TDK (IA281 . and (IA291 one obtains the expression for V / (-|+,-) as given in Eq. d20l >. Similarly, after rewriting 
Eq. dAUb for > 1 as 



■=■4 



Ai>(S,0>l,A) = (A_-A + )0 J e-®^- f ds fjarccosfa -1 / 2 ) - V^~l") tr 3 * */ 2 , 
one can integrate by parts, which yields 



A^(S,0>l,A) = (A_-A+)0 rf -V / ds—^= 

J\ \/s— 1 



1 0S 2 



s 2 
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(A31) 
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Using Eq. ( !A27b and the relations [which follow from taking successive derivatives —d/d{a 2 ) of Eq. ( IA27M 



\ 1 -ah VX -a 2 

as , -.e = e 



-a l s V JL /, r, 2\ -a 1 



ds- 



one ends up with 

A#(S,0> 1,A) = 7r(A_ -A + )0^ 2 e- erfc(S v / 072) 
Together with the expression for the homogeneous case [see Sec.imiand Ref.uD, 



d-2-@ 



% -) (® > 1 , A -> 0) = 2^A ± fl - z e 



(A32) 



(A33) 



(A34) 



one obtains the expression for fi)(_i+ -) given in Eq. d20b . 



Appendix B: Derjaguin approximation for a single chemical lane 

Based on the assumption of additivity which underlies the Derjaguin approximation one can use the results presented in Sec.lIVI 
for a chemical step in order to study a chemical lane. The chemical lane configuration can be regarded as the superposition of 
two chemical steps, one (A) being a chemical step located at x — — L with {a\ai) BC, and the other one (B) being a chemical 
step located at x — L with (ai\a) BC. This superposition overcounts a contribution corresponding to a homogeneous substrate 
with (of) BC which must be subtracted [see Eq. (fTJtl: 
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K, 



B 



(a,b) + K (a e ,b) 



(.at) 



( fl > Jfj^l+JOL 



-L 



K, 



(a h b) 



K (a,b) ~ K {a t ,b) 



Kr 



^( fl |^,6)(S+A,0,A) 



(Bl) 



(B2) 



and 



4 aj (A,S,0,A) 



A", 



{a,b) + K (a e ,b) 



K (a t ,b) - K (a,b) 



Vteh*)(S-A,0,A). 



(B3) 



Since within the DA y/( a ,i a ^\ = W(a\a,,b)^ Eqs. ( IBU - (IB3b and Eq. ( f27l > lead directly to Eq. ( f29b . The procedure for calculating 
the critical Casimir potential is analogous to the one discussed here for the force and leads to Eq. d30b . 



Appendix C: Derjaguin approximation for periodic chemical patterns 

In order to obtain the scaling function for the critical Casimir force and the potential of a sphere close to a periodic chemical 
pattern one can follow a procedure analogous to the one presented in Appendix|B] Indeed, in order to form a lane £' with (a\) 
BC on an otherwise homogeneous portion of a substrate with (02) BC and lateral extension P, one can proceed as follows: 

(A): superimpose onto the substrate the single chemical lane t studied in Sec.JVJ with ag =a\, a = ci2, suitably positioned in 
space such that it coincides with the lane £' to be formed. 

(B): subtract the contribution of a homogeneous substrate with (02) BC, which is overcounted in the previous superposition. 
After this subtraction, the contribution to the force resulting from that part - marked by (?) in Eq. ( IC1I ) - of the original 
substrate which is not affected by the formation of the extra lane is unchanged. 
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(?) 



M 



-nsmsp 



(A):0- 



(«2) 



(»l) 



(«2) 



--J- A +T 

(«2) 



(CI) 



lane £' 



The contribution AF to the critical Casimir force experienced by a colloid close to such a substrate and due to the addition of 
the lane is characterized by the scaling function [see Eq. d27l il 



K (a 2 ,b)- K (a,,b) 



AK(X > n > E-E',®,A^0)=Kt(n^E-E\®)-K (a2tb) 



V(a 1 |« 1 »)(s-s'+n§,®,A-+o)-vf(« 1 | fl| b)(s-s'-n§ ) 0,A->o) 



(C2) 



where we have used the relation (Li/2)/\/RD = TlX/2 and have introduced S' = X'/^/RD, with X' as the position of the center 
of the added lane (! . The force resulting from a periodic pattern can now be obtained by starting out with a homogeneous 
substrate with (02) BC and by iterating the procedure discussed above which adds progressively displaced lanes at positions 
X' = nP, i.e., S' = nil, with n e Z. The resulting force is characterized by the scaling function 



+°° 



K p (X,n,z,e,A^o)=K {a2 j }) + £ A^(A,n,s-«n,@,A^o) 



(C3) 



which, together with Eq. d33l ). yields immediately Eq. ( l35l l for i// p . 

For A = or X = 1 one recovers from Eq. (l35l l the homogeneous cases with (02, b) BC or (a\,b) BC, respectively. Obviously, 
for X = 0, the sum in Eq. d35l l vanishes, and one is left with y/ p (A = O,n,E,0,A-> 0) = 1, corresponding to (02, b) BC. On the 
other hand for A = 1, the sum in Eq. d35b can be easily evaluated [see Eq. ( TTol ) for |S| — > °°]: 



M.N 



1 jP E {^(« 1 |« 2 ,2 ) )(S+n(n + i) ) 0,A)-VA (fll | a2)6 )(S+n(n-i),0,A)} 



= -M 



= JJ£ {^(a 1 | a2! fe)(S + n(Ar+i),0,A)-VA (ai | a2ife) (S+n(-M-i),0,A)}=-2, (C4) 



M,N- 



where we have used the fact that Y(a 1 \a 2 ,b){'^ = ±°°,®,A) = =pl. Accordingly, y/ p (A = 1,II,S,0,A — > 0) = — 1, which corre- 
sponds to the homogeneous case with (a\,b) BC. 

In the limit II — » (i.e., for very fine patterns compared with ^/RD), the sum in Eq. d35l l turns into an integral: 

E {^( ai |a 2l 6)(S + n(n + |),0,A)- V A (aiK>6) (S+n(n-|) ) 0,A)} 

n.=— 00 

= J _J*\ X ^V / i>ik!,2>)( E + T l> e) > A ) 

= ^ { V( ai |W+°°' ' A ) - ^ (fll | a2)6) (-oo,0,A)} = -2A, (C5) 

and finally one finds Eq. (l36l l. 

For completeness, we provide the corresponding expression for the scaling function of the critical Casimir potential C0p within 
the DA: 

C0p(A,n,S,0,A^O) = l+ £ {a) Wfl2i6) (S+n(n+|),0,A^O)-fi) Na2>fe) (S + n(n-|),0,A^O)}. (C6) 

n=— °o 

In the limit II — > 0, Op reduces to 



a)p(A,n^O,S,0,A^O) = 1-2A. 



(C7) 
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Accordingly, within the DA and in the limit n — >• the critical Casimir potential is the average of the ones corresponding to the 
two boundary conditions, weighted with the corresponding relative stripe width: 



v p (X,n -► o,s,e, a -> o) = x$ (aub) (&,A -> o) + (1 - A)0 (fl2i6) (e, a -». o). 



(C8) 



Appendix D: Cylinder close to a patterned substrate 



1 . Derjaguin approximation for a homogeneous substrate 



Similarly to the case of a sphere discussed before, the crit- 
ical Casimir force Ff y ,< per unit length acting on a (three- 
dimensional) cylinder of radius R with (b) BC close to and 
parallel to a substrate with (a) BC at a surface-to-surface dis- 
tance D can be expressed in terms of a universal scaling func- 
tion K Q y l : 



F Z)(D,R,T)=k B T 



fll/2 



D d-l/2^(a,b) 



** (0,A), (Dl) 



with © = sign(f ) D/E,± and A = D/R as before. Equation (IDU 
describes a force divided by a length and per Z/'~ 3 which for 
d = 4 corresponds to considering F? y b -> per length of its axis 
and per length of the extra translationally invariant direction 
of a hypercylinder [compare Eq. ©]. The geometric prefac- 
tor in Eq. dDU . however, differs from the one for the sphere 
[Eq. (0] because it is chosen such that within the DA (A — > 0) 
the scaling function K? y , > attains a nonzero and finite limit, as 
discussed before. The DA can be implemented along the lines 
of Sec. MI Al for the sphere. Here the surface of the cylindrical 
colloid is decomposed into pairs of infinitely narrow stripes of 
combined area dS = 2Mdp , positioned parallel to the substrate 
at a distance L(p) from it [Eq. ^}] and each at a distance p 
from the symmetry plane of the configuration. M is the length 
of the cylinder and drops out from F ? y , > which follows analo- 



gously from Eqs. ([8]) and (O: 



-,cyl 



u.b) 



Fg l fi) (D,R,T)/k B T~2j o dp[L(p)}- a k (aM (sign(t)L(p)/^) 

(D2) 
where L(p) is given in Eq. (0. Finally, in the limit A — > we 
obtain 



*?>.A- 



0) = V2 / da (a- 1)~2 a ~ d k (a . b) {®a) 

i 



At the 

\a,b)^ 



^ M (0,0) = 



bulk critical point 

2n\T{d - i) 



K?' (0,0) = [3^/(4V5)]A (a 





^cyl 



^ ) (0,0) = [5w/(8V5)]A ( , 



(a,6) 





(D3) 


= one 


finds 


>/r(d)]A (0]6 ) so 


that 


1 .66 x A( fl ^ for d = 


= 3 and 


1.38 x A( ai j for d = 


= 4. 



2. Derjaguin approximation for a chemical step 

Here, we assume that the axis of the cylinder is parallel to 
the chemical step, i.e., perpendicular to the x direction [Fig.[T|, 
as well as parallel to the substrate. The projection of the posi- 
tion of the axis of the cylinder with respect to the x direction 



I 

is denoted by X, so that at X = the cylinder is positioned 
directly above the chemical step [Fig. [T]- Accordingly, the 
problem is effectively two-dimensional and the correspond- 
ing DA can be performed much easier than in Appendix [A] 
Following an approach analogous to the one adopted for the 
sphere in Sec. [IV] and in Appendix lAl we rewrite the normal 
critical Casimir force per unit length acting on the cylinder as 
in Eq. dA3l ): 

F s cy] (X,D,R,T)=F{ yl >h) (D,R,T)+AF cy \x,D,R,T). (D4) 

Within the DA we find for A — > [compare Eq. ( IA6H 

AF cy \x,D,RJ)=k B T-^ I - m AK cyl (Z,®,A^0), (D5) 
where [compare Eq. ( IA7M 
A^ cyl (E,©,A^0) = \/2 / da(a-lHcT d M(0a). 

Jl+S 2 /2 

(D6) 
Using Eq. dD6b and Eq. dD3b we find for the whole range of 



,, c y' 



which is defined 



values of S the scaling function w, 

completely analogous to Eq. ( [T3l [compare Eq. (TTSll: 



,,cyl 



V^| fl>1 *)(S^0 J ® > A-^0) = Tl 



± 



V2 Jr+z2/2 da («-!)"' a~ d Ak(& 



a 



<I^m(©>A^O)-^ ) (0,A^O)' 



(D7) 



(a<,bV 



3. Derjaguin approximation for a periodic chemical pattern 

The derivation of the scaling function for the critical 
Casimir force acting on the cylinder close to and aligned with 
a periodic chemical pattern as studied in Sec. lVHI is analogous 
to the one for the sphere described in Appendix [C] The final 
formula for \f/p Y is the same as in Eq. ( 1351 ) with yff ai \a 2 ,b) re " 
placed by y^ , . . given by Eq. ( 1D7| |. This renders the critical 
Casimir force per unit length 



F' yl (L u P,X,D,R,T) 



hTj^-^K; y \X,Tl : Z,® : A) 



where K^ y is defined as in Eq. d33l l with Kr a b -\ and K^ 



(D8) 

a 2 ,b) 

replaced by K~/ a ' b) and ]££ b y respectively, which are given 

by Eq. dP3t , and with y/ p replaced by y/p y . The corresponding 
results are shown in Fig.[T2l 



jrcyi 



^cyl 
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